Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 15 December 201 1 (MN MbX style file v2.2) 



Surface Brightness and Intrinsic Luminosity of Ellipticals 



Barun Kumar Dhar^*, Liliya L.R. Williams^ 

^School of Physics and Astronomy, University of Minnesota, 116 Church Street SE, Minneapolis, MN 55455 USA 



Accepted 1 2th December, 20 11 ; in original form 2nd June, 2011. 



ABSTRACT 

We show that the surface brightness profiles of elliptical galaxies can be parametrized using 
a linear superposition of two or three components, each of which is described by functions 
developed in Dhar & Williams as the 2D projections of a 3D Einasto density profile. For 
a sample of 23 ellipticals in and around the Virgo Cluster with total absolute V-magnitude 
— 24<A/yT< — 15, our multi-component models span a dynamic range up to 10® in surface 
brightness and up to 10^ in radius down to the resolution limit of the HST, have a median rms 
of 0.032 mag arcsec~^ consistent with the rms of 0.03 from random errors of the data, and 
are statistically justified at >3cr. Our models indicate that i) the central component is more 
concentrated than the outer component; and ii) the central component of massive shallow-cusp 
('core') galaxies is much more luminous, extended and concentrated than that of steep-cusp 
('cuspy') galaxies, with their near exponential central profiles indicating disk-like systems, 
whose existence, must be verified spectroscopically. 

Galaxy structure can thus be modelled extremely well with a central mass excess for all 
galaxies. This is not necessarily contrary to the notion of a mass deficit in 'core' galaxies, 
since mass ejection due to core-scouring by a supermassive black hole (SMBH) binary could 
have affected the shape of the central components. However, we show that the existence, 
amount, radial extent and sign of such deficits disagree substantially in the literature, both for 
a given galaxy and on an average over a sample. We discuss possible implications and suggest 
that SMBH binaries are unlikely to be the sole mechanism for producing the large 'cores' of 
massive galaxies. 

Using results from the SAURON survey we deduce that under certain conditions of sym- 
metry, inclination angles and degree of triaxiality, the intrinsic (3D) density of light can be 
well described with a multi-component Einasto model for both steep- as well as shallow- 
cusp galaxies. This indicates an universality in the functional form describing the 3D density 
distribution of light in galaxies and dark matter in ACDM N-body haloes. Finally, Plane- 
tary Nebulae and strong lensing observations, and the Einasto index n of ACDM dark matter 
haloes indicate that our result - the outer component of the surface brightness profiles of mas- 
sive galaxies has 5 < n < 8 - could imply i) a common feature of collisionless systems; and 
ii) that galaxies with such n for their outer component are dark matter dominated. 

Key words: galaxies: structure - galaxies; photometry - galaxies: fundamental parameters - 
galaxies: haloes - dark matter - gravitational lensing 



1 Introduction 



iGerard de VaucouleursI ( Il948h showed that a remarkably simple 
parametrization of the radial surface brightness profile of galaxies 
exists for a wide range of ellipticals. He proposed a two parameter 
fitting function of the form of equation i ll. It with m=A defining 
the shape of the distribution and q—7.67 ensuring that the effective 
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or half-light radius. Re, contains half the total projected light. 



E(7?) = T.R^ exp i -g 



R 

Re 



- 1 



(1.1) 



where, E(7?) is the 2D surface brightness at a plane of sky projected 
radial distance R, T,r^='E{Re), q—q{m) and Es(0)=E_R^e''. 

Observin g that our Galaxy is made up of multiple subsystems 
lEinastdd 19651) proposed a modification, equation l ll.2t . of the 2D 
de Vaucouleurs law to model the intrinsic (3D) baryonic mass den- 
sity p(r) of each subsystem by allowing the shape parameter n to 
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be a free parameter, and b—b{n). 

p{r) = ps exp < -6 



- 1 - 1 



(1.2) 



where, ps is the density(3 D) at a scale ra dius and p(0)=pBe'' . 

Around the same time lSersi3 l ll968h observed that m in equa- 
tion l ll.lt . characterized by m = 4 in de Vaucouleurs law, was not 
the same for all galaxies and rendering it a free-parameter provided 
much better fits to the surface brightness profiles (hereafter, SB). 
Equation Jl.U is the current standard paradigm for describing the 
global structural properties of galaxies. 

Over the past forty years, the pioneering works of Caon, Ca- 
paccioli, Einasto, Ferrarese, Graham, Kormendy, Lauer and their 
collaborators have shown that no single commonly used three pa- 
rameter fitting function could model the SB over the entire dynamic 
radial range. They showed that the SB profiles of all galaxies re- 
veal an inherently multi-component structure such that the outer re- 
gions can be described with a Sersic or a Sersic-Fexponential mod- 
els while the central regions can be described with power-laws. 

Additionally, for some nearby spirals and the giant elliptical 
M87 (NGC4486), Einasto and collaborators (Rummel, Haud and 
Tenjes) have shown that if spectroscopic and kinematic constraints 
are used in addition to the SB data, then the intrinsic 3D mass den- 
sity including their central regions can be described with a multi- 
component Einasto model. 



1.1 Motivation for tliis work 

While the power-law-l-Sersic models are widely accepted as pro- 
viding an accurate description of the SB profiles, the fit residuals 
of these models are often larger than measurement errors. Further, 
these functions do not have the flexibility to model deviations from 
power-laws within the central regions. Since these fits are often 
used to draw inferences on galaxy structure, formation and evo- 
lution, it is important to address two crucial issues when modelling 
galaxy structure: 

(i) Models of the SB profiles must be consistent with measure- 
ment errors over the entire available radial range; and 

(ii) Model fitting functions must allow one to easily infer the 
intrinsic 3D luminosity density from the 2D SB profiles. 

Our goal in this p aper is to address the ab ove issues with a 
new function derived in lDhar& Willi"amsl ( l2010l) (DWIO) to model 
the projected surface mass density profile of Einasto-like systems. 
We show that a multi-component parametrization with this function 
(hereafter, the DW-function) provides excellent fits, consistent with 
measurement errors, to the SB of ellipticals over a large dynamic 
radial range down to the HST resolution limit. 

This suggests that the 3D density distribution of light in galax- 
ies can be described with a multi-component Einasto model. 

Such an interpretation is similar to that of Einasto and collab- 
orators, except that a) our work extends to a much larger sample of 
shallow and steep cusp ellipticals; and b) we deduce the intrinsic 
3D Einasto model properties by fitting the 2D SB using a multi- 
component DW-function. 

For the rest of this paper, we draw a distinction between com- 
ponents deduced from fitting and physically distinct kinematic sys- 
tems or stellar populations. We shall refer to a component as a sin- 
gle spherically symmetric closed form fitting function described by 
at most three parameters - a scale length, r^, a shape parameter, a, 
and some normalization, Eq or Eq. Our multi-component fits con- 
sist of the minimum number of such functions (i.e. DW-functions) 



required to model the entire dynamic radial range down to the res- 
olution limit of the HST. The minimum number depends on the 
quality of data, the available degrees of freedom, and the amplitude 
of systematic patterns in the residuals. 

One should hence use caution in interpreting these compo- 
nents as kinematically distinct systems or stellar populations since 
the physical properties of the components, and even their presence 
will depend on the choice of parametrization. It is however possible 
that some of our fitted components do coincide with kinematically 
identified systems or stellar populations, or combinations thereof, 
which shall then facilitate interesting interpretations. 

Also note that neither fit components nor physically distinct 
systems have to correspond to structure in the total gravitational po- 
tential. Only where pbaryons':>PDM , or pbaryonsCcpoM fcatures 
observed in the SB profiles will trace the total mass density, and 
hence will provide information about the gravitational potential. 

Further, since no galaxy has a truly flat density core, we refrain 
from using the terms 'core-' and 'cusp-' galaxies and instead refer 
to them as shallow-cusp and steep-cusp galaxies, respectively. 

1.2 Brief outline and summary of results 

In section [2] we provide a brief history and basis for believing that 
galaxies have a multi-component structure. This is followed by an 
overview, in section [3] of the most flexible fitting functions tried 
to date, and the motivation for this work. The data we use are de- 
scribed in section |4] and important features of our 2D parametriza- 
tion and the fitting procedure are discussed in section |5] Section 
[6] discusses results of our 2D fits to a sample of 23 ellipticals in 
and around the Virgo cluster spanning absolute V-magnitudes in 
the range -24<Mvt<-15. Our models produce consistently low 
residuals, over a large dynamic radial range (~10^) and in congru- 
ence with the measurement errors. 

Section|7]contains a comparison of fits with other parametriza- 
tions from the literature. It highlights the main result of this work: 
of all the functions tried to date to model the 2D structure of galax- 
ies, the best fitting model is the multi-component DW-function. This 
suggests that the intrinsic 3D density structure of galaxies is best 
represented with a multi-component Einasto model. 

In section [8] we discuss the properties of the components 
deduced from the multi-component fits with the DW-function, 
and in section [8!4l we present three cases where the central DW- 
component of our multi-component fits coincide with spectroscop- 
ically identified systems. 

We next explore the conditions under which our 2D models 
can be used to infer that the intrinsic 3D luminosity density dis- 
tribution can be described with a multi-component Einasto model. 
In section [9] we therefore first discuss the key issues pertaining to 
non-uniqueness of deprojection, and then in section[TO]we present 
our models of the 3D Einasto luminosity density profiles for 14 
galaxies whose deprojections are less likely to suffer from non- 
uniqueness. Since these 14 galaxies span a wide luminosity range 
and belong to both the steep- and shallow-cusp classes, we con- 
clude that the intrinsic structure of all galaxies could be described 
with multi-component Einasto models as well. 

This is followed by section [TT] where we discuss the notion 
of 'mass deficit' in massive ellipticals. Following a brief review of 
the phenomena of mass ejection by binary SMBHs and results from 
N-body simulations, we show that estimates of the amount, radial 
extent and sign of the 'observed' deficit in real galaxies, have large 
disagreements in the literature - both for a given galaxy as well as 
on an average for a given sample. 
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In this paper we have shown that the structure of galaxies can 
be modelled extremely well with a central mass excess for all galax- 
ies. We suggest that, while mass ejection due to core-scouring by 
binary SMBHs could have shaped the central and the intermediate 
DW-components, this phonemena by itself is unlikely to be the sole 
cause for the observed shallow-cusps; other formation processes 
are likely to have also contributed to the existence of this feature. 

In section[T2]we show that an Einasto shape parameter of 5 < 
n < 8 is found in two different systems - (i) the outer baryonic DW- 
component of the massive ellipticals in our sample, and (ii) in fits 
with an Einasto profile to the collisionless ACDM N-body haloes. 
This range of n may thus be a common feature of collisionless 
systems. Finally, by comparing results from Planetary nebulae and 
strong lensing studies we show that it is possible that galaxies with 
5 < n < 8 for the outer component of their SB profiles are likely to 
be more dark matter dominated. 

Section [T3] contains a general discussion and a summary of 
conclusions drawn from the work presented in this paper. 



2 Multi-component structure of galaxies 
2.1 Structure of baryonic (stellar) distribution 

Hubble's (1926) classification of galaxies as ellipticals and spi- 
rals based on the identification of bulges and disks was probably 
the first systemat ic characterization of the multi-component struc- 
ture of galaxies. Ide Vaucouleurs showed that the outer 
regions of spirals and lenticulars can be modelled as an expo- 
nential (also see Patterson (1943), de Vaucouleurs (1955, 1959b) 



cusp ('core') el liptical M87 (N GC4486), Einasto and collabora- 
tors (following lEinastd ( Il965h ) demonstrated that the 3D bary- 
onic mass density can be described as a multi-component sys- 
tem of nested Einasto profiles by parametrizing each kinematically 
and photometricall y identified system o r stella r population using 



with the central bulge following an exp( — r ' ) profile, similar 
to elli pticals (de Vaucouleurs 1948,1958). iBertola & Capaccioli I 
modelled the giant elliptical M87 as a linear superposition 
of a de Vaucouleurs and exponential profiles. The multi-component 
structure of lenticulars as a superposition of an exponentia l and a 
exp (— r^^"^) compone nt was also investigated in lFreemMil l ll970l) 
and lKormendy] l ll977h . 

However, it was also observed that not all ellipticals and 
bulges of lenticulars followed a pure de Vaucouleurs profile. Even 
then Sersic's (1968) generalized form, equation dl.lt . did not find 
much appl ication until the com prehensive works of ICaon et al] 
( 1 19931) and lGraham eTal] ( Il996h . which revealed how an incorrect 
parametrization - a de Vaucouleurs profile instead of a Sersic pro- 
file - could lead to misleading con clusions (also see iDavies et al.l 
( Il998l) and lYoung & Curi3 ( 1 19941) ). Since then the Sersic pro- 
file has become the norm for describing t he SB of galaxies ove r 
large dynamic rad ial ra nges. Additionally, Ferrarese et alj l ll994h . 
iLauer et alj ( 1 19951) and lGraham et al. I ( l2003h also observed that a 
single Sersic profile cannot fit the entire dynamic radial range, es- 
pecially the central regions down to the HST resolution. The failure 
of a single Sersic profile to fit any of the galaxies in their sample 
confirmed that multiple fit components are necessary to adequately 
model the structure of galaxies. 

In the absence of a sound theoretical model of galaxy forma- 
tion and evolution it is important to understand under what condi- 
tions can analysis through fitting functions identify the true intrinsic 
components of a galaxy. This is especially important since a correct 
(or incorrect) identification of such components can have a serious 
bearing on our understanding of galaxy structure, formation and 
evolution. It is therefore important to ask, how many components 
in a galaxy can be realistically identified and how? 

In a series of papers on nearby spirals and the giant shallow- 



equati o n l|1.2|> (see lEinasto & Rummel I |l970l). iHaud & Einasto 



1989 ), Tenies. Einasto & Haud (1199 II) . iTenies. Haud & Einasto 
19941) . iTenies Haud & Einasto \ ( Il99i) ). They observed that, de- 



pending on the galaxy and quality of data, the 3D mass density can 
be described as a certain combination of superposed systems - a 
nucleus, a core, a bulge, a disk, a halo, a flat system - all well fit 
with equation M.2\ and a massive corona of baryonic and dark mat- 
ter described using an isothermal density profile. It is to be noted 
that in their construction each Einasto-component had an unique 
shape parameter n and scale-length Ts identifying physically dis- 
tinct systems as opposed to apparent components generated by the 
use of ( ll.2t merely as a fitting function. To our knowledge, this was 
the first systematic modelling of the intrinsic 3D structure of the 
components in galaxies through a prior identification of systems. 

In the pre-HST era, while the detection of bulges and disks 
revealed a multi-component structure at larger radial scales, the 
seeing effects of the ground based observations did not allow one 
to adequately resolve the central regions. Some galaxies possibly 
showed central flattening of the density profiles, but because of lim- 
ited resolution these claims remained inconclusive. 

The dawn of the HST and high resolution, 0.02 arcsec 
pixel~ ^, imaging improved the situation dramaticallv. ICrane et al.l 
( Il993l) . using the HST Faint Object Camera (FOG) showed that 
the central regions of galaxies do not have flat density cores, but 
instead exhibit a wide range of inner power-law slopes. The low 
S/N of their observations however prevented them from obtaining 
psf-deconvolved images that would allow for detailed modelling. 

Imaging wi th the Wid e -Field Planetary Camera I ( WFPC - 
I) on the HST IjaffeetalJ jl994h. [van den Bosch et al.l d 19941) . 
iFerrarese et alj jl994h and iLauer et al. ( 19951) generated psf- de- 
convolved images of the central regions of ellipticals. They found 
that the SB of the central regions convincingly revealed a wide 
range of inner slopes inside a characteristic break, or transition 
radius, with shallow slopes and fainter central magnitudes for the 
largest and most luminous galaxies, and steeper slopes and brighter 
central magnitudes for the smaller and less luminous galaxies thus 
implying an inherently diverse multi-component structure. This led 
to the development of a num b er of power- l aw bas e d parametriza- 
tions jFerrarese et all jl994h . iLauer et al] jl995l) . iGraham et all 
(I200I)) to model the central and outer regions of galaxies. 

Thus analysis of the SB profiles through parametric fitting 
functions also revealed a multi-component structure for galaxies, 
indicating the presence of two or three regions; a double power-law 
domain around a transition radius, and a main body for the galaxy 
usually well described with a Sersic profile. This is similar to the 
work done in modelling the 3D mass density by Einasto and co- 
workers except that there is no reason to assume that components 
generated through the use of fitting functions should correspond to 
kinematically identified components or distinct stellar populations. 
Further, there is no convincing reason as to why the form of the 
fitting functions should be different in different regions. 

2.2 Structure of ACDM haloes 

iNavarro et al.l ( |2004|) (N-l-04) showed that 3-parameter fitting func- 
tions, especially those with a power-law logarithmic slope like the 
Einasto profile, are able to describe the 3D mass density distri- 
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bution of spherically averaged ACDM N-body haloes better than 
any of the 2-pa rameter fittin g funct ions tried to date. Subsequen t 

l |2006|) . 



simulation s by 

jgOOgj , 



Gaoetal 



'Merritt et al.' (2006) (M+06), Prada et al 
HayashL& White (2008), Stadel et al 
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Navarro et aiT i lio id) have confirmed the need for a 3-parameter 
function and for over 30 such dwarf, galaxy and cluster size N- 
body haloes the Einasto profile seems to be the best performing 
function in comparison to other 3-parameter fitting functions. The 
Einasto index deduced from these simulations are typically in the 
range 5<n<8. IStadel et alj ( |2009|) proposed a 2-parameter func- 
tion that provides fractionally better fit in terms of rms than the 
Einasto profile for two haloes they simulated. However, even for 
these cases the Einasto profile has comparable residuals. 

Note that current resolution of N-body simulations does not al- 
low one to probe the very central regions of galaxies where baryons 
reveal a multi-component structure. One may thus conclude that 
within the resolved and converged domain of N-body simulations, 
the dark matter distribution can be described as a single component 
system. 

2,3 Intrinsic and Projected structure 

Noting the similarity b etween the fu n ctiona l form of the Einasto 
and the Sersic profiles, iMerritt et al.l l l2005h obtained spherically 
averaged non-parametric estimates of the 3D intrinsic and 2D pro- 
jected mass densities of the N-body haloes described in N-l-04. They 
found that the same fitting function, the Sersic profile, which is used 
to describe the SB of ellipticals also describes the surface densities 
of ACDM N-body haloes, whose intrinsic 3D density is best de- 
scri bed by a function of simil ar form - the Einasto profile. 

iDhar & Williamsl J20T0I) show that while it is possible to find 
limited radial ranges over which a Sersic profile can approximate a 
projected Einasto profile, over a large radial range a Sersic profile is 
not a good representation of a projected Einasto profile, and using 
such fits can lead to a misinterpretation of the best fit parameters. 
DWIO point out that the fitted Sersic profile parameters depend 
strongly on the radial range of a projected Einasto profile. In other 
words, fits with a Sersic profile to a single projected Einasto profile 
implies the existence of a variable Sersic shape parameter. DWIO 
provide an accurate analytical approximation for the 2D projection 
of an Einasto profile in terms of the 3D Einasto profile parameters, 
thereby allowing one to explore the intrinsic properties of systems 
believed to be Einasto-like from 2D observations of those systems. 



3 An overview of fitting functions 
3.1 Combination of power-law functions 

With the ability to resolve the central regions of galaxies with 
the HST, it was found that a single 3-parameter fitting function, 
like the Sersic profile, is not able to model the SB profile over a 
large dynamic radial range down to the HST limit; less so with 2- 
paramet er functions, for exam ple, the de Vaucouleurs and the Jaffe 
profiles. [Ferrarese et alj 1 1994h thus prop osed the 4-parameter dou- 
ble power-law while Lauer et al] i ll 9951) proposed the 5-parameter 
Nuker profile - which is a modified double power-law with an addi- 
tional parameter to control the sharpness of transition between the 

pow er-laws. 

iLauer et alj ( Il995h pointed out that the more flexible 5- 
parameter Nuker profile is designed to model only the small cen- 
tral regions, ~ 10-20 arcsec, and not the rest of the large galaxy 



structure, which could be well fit with a 3-parameter Sersic profile. 
This is because power-laws assign a fixed logarithmic slope to the 
density distribution while the light of galaxies at large R exhibit a 
variable slope. However, even with 5 parameters modelling a small 
radial range it can be seen that the central-most regions show resid- 
uals larger than measurement errors. This implies that a total of 
8 parameters are required for a near complete description exclud- 
ing the very central inner regions. Moreover, in the central regions 
many galaxies exhibit a sharp change in slope, i.e. a transition ra- 
dius. Even though the presence of the transition is unambiguous, 
the domain over which one should fit the Nuker profile to obtain 
robu st parameters is not a lways obvious. 

lOraham et alj l l2003h showed that the best-fitting Nuker pa- 
rameters are extremely sensitive to the chosen domain of fit. Hence, 
observing that the central regions can be modelled as power-laws 
while the outer regions require a Sersic profile, they proposed the 
6-parameter Core-Sersic profile which has an inner power-law cou- 
pled to an outer Sersic pr ofile along with a pa rameter to control 
the sharpness of transition. iTruiillo et alj l l2004h suggested a mod- 
ification of the Core-Sersic profile by allowing an infinitely sharp 
transition along with a step-function that reduces the Core-Sersic 
profile to a 5-parameter function, but as expected produces an un- 
physically sharp break in the profile quite unlike the much smoother 
transition observed in galaxies. Even then, the large fit residuals in 
the central region continue to exist. 

iFerrarese et alj 1 I2OO6I) found that the light excess in the cen- 
tral regions (often referred to as 'nuclei') of the steep-cusp galaxies 
in their sample could be best described with a 3-parameter King 
model. Such 'nuclei' are not as rare a s was originally believed; 
but may not always be very prominen t. ICote et alj 1I2OO6I) showed 
that in the ACS Virgo Cluster survey l lCote et alj|2004l) . hereafter 
ACSVCS, 66-82 per cent of the galaxies have such a central fea- 
ture. It thus appears that one needs a King model for the nuclear 
region, and either a Sersic, a Core-Sersic or a Nuker-l-Sersic profile 
to model the rest of the galaxy, i.e., a total of 6-1 1 parameters. 

While the overall rms of the fits with 6-11 parameters appear 
small, the residuals in many radial sections remain large, > 0. 1 
and sometimes > 0.2 mag arcsec~^; considerably larger than the 
0.01-0.05 mag arcsec^^ uncertainty of HST qua lity data. For il - 
lustration we refer the reader to the fits in Ferrares e et al. I j2006l) 
using a combination of Sersic, Core-Sersic and King models where 
large, and sometimes divergent, residuals can be found in some 
regions. A comparison of fits with Core-Sersic, King and Nuker 
profiles in the cen tral regions of some galaxies are also shown in 
iLauer et al] ( |2007|) . clearly revealing the lack of a good fit in the 
central-most regions. More detailed discussion is provided in ij?] 

Further, none of the above power-law based functions are de- 
fined in the limit R ^ 0. While this may seem to be of academic 
interest, since the density at i? = does not contribute much to the 
light enclosed, it can introduce uncertainties in interpreting depro- 
jections needed to extract 3D intrinsic profiles, and also for fitting 
procedures using psf-convolved models, where one needs to spec- 
ify a finite value for the SB at i? = for the convolution. 

3.2 Multi-component modified exponentials 

All of the above power-law and Core-Sersic parametrizations are 
well-guessed but ad-hoc empirical fitting functions in 2D, in the 
sense that they are not a result of well established theoretical mod- 
els of galaxy formation. However, despite the existence of residuals 
larger than measurement errors, the parameters of the fits are used 
to draw inferences on galaxy structure and evolution, which are in- 
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trinsically 3D phenomena. This could have been meaningful if the 
2D models were deduced from an underlying physically motivated 
3D distribution. In order to draw such inferences from functions 
that are merely fitting functions, one needs, at the very least, to 
have residual profiles consistently comparable to measurement er- 
rors over a large dynamic ra dial range, and not just a low rms. 

iKormendv et al.l ( l2009h . hereafter KFCB09, addresses this is- 
sue partly by fitting a single modified exponential function, the Ser- 
sic profile, over a rigorously tested range where the function pro- 
duces residuals comparable to measurement errors. However, this 
range is chosen to ensure that the Sersic profile produces a good fit 
and specifically excludes the entire domain within the break radius. 
Since the Sersic profile is an ad-hoc fitting function, parameters de- 
duced from fits within a limited radial range can lead to misleading 
interpretations of the physics involved and in estimating properties 
of the region inside the break radius. 

Since the Sersic profile produces good fits over a large radial 
range and lenticulars have been modelled with a Sersic -I- exponen- 
tial profile, a plausible alternative is to use a double-Sersic profile 
to model the entire galaxy. This approach had not received much 
attenti on until recently when iGonzalez. Zabludoff & Zaritskv I 
l l2003h showed that a double-Sersic profile provides much better 
fits to twelve Brightest Cluster Galaxies (BCG). An example of 
such a fit for the BCG in Abell 2984 is shown in their figure 2. 
ICoteetal.l ( l2007l) also arrive at a similar conclusion for some galax- 
ies in the ACSVCS and Fornax cluster surveys. A comprehensive 
study of fits with a double-Sersic profile is provided in Hopkins et 
al. (2009a,b) where they obtain low residuals over a large dynamic 
radial range. While they do not provide a residual profile, the rms 
of their fits are often larger than the measurement errors. 

From a mathematical stand point, the Sersic profile presents 
another difficulty: its deprojection, that can provide insights 
in to the 3D structure of the galaxy, is not very well an- 
alytically tracta b le. A symptotic limits of deprojection can be 
fo und in Ciottil 1 1991) and appr oximate expressions are given 
in IPrugniel & Simiei] i 1997h and iLima Neto. Gerbal & Marguezl 
( 1 19991) (hereafter. PS97 and LGM99). However, the PS97 and 
LGM99 approximations are not accurate at small R^IQ'^Re, and 
the 3D densi ty becomes undefined as r— >0 for Sersic indices m>l 
( ICiottill99lh. while all galaxies observed to date have m>l. 



Baes & Gentile (BGll) provide an exact analytical 



expression for the deprojection of Sersic profiles for all m in terms 
of the Fox H function, which is extremely difficult, if not impossi- 
ble to compute, even numerically. However, for rational values of 
m they show that the deprojection can be expressed in terms of the 
Meijer G function. Rational m requirement is not too stringent be- 
cause for practical purposes any m can be well approximated by a 
rational number. The singularity in the deprojected central density 
for m > 1 is, however, inherent to the form of the Sersic function. 
We also note that the deprojections discussed above are assumed to 
extend to infinite 3D radius. However, just because a deprojection 
is analytically difficult, does not, by any means, suggest that the 
true 2D light distribution of galaxies is not described by a Sersic 
profile, and the 3D distribution is not a deprojected Sersic profile. 
What might suggest that the Sersic profile is not an optimal function 
over very large radial ranges is that even a 2-component model, i.e. 
a double Sersic (as in Hopkins et al. (2009a,b)), often yield residu- 
als larger than measurement errors. 

Motivated by the finding that the Einasto profile provides bet- 
ter fits to pure dark matter high resoluti on N-body simulat i ons of 
dwarf, galaxy and cluster sized haloes, iDhar & Williams! | |2010|) 
presented, for the first time an extremely accurate - fractional de- 



viations of ~ 10~* to 10~^ - analytical approximation to the sur- 
face mass density of a 3D Einasto profile. This function is valid 
for n > 0.2 (see section (|5} below), and is expressed in terms of 
the 3D Einasto profile parameters, ps, Ts, and n. Given the issues 
described above with the existing forms of the fitting function, in 
this paper we explore the quality of fits to the surface brightness 
profiles of ellipticals with a multi-component DW-function, which 
has the interesting property that the intrinsic 3D luminosity density 
is a multi-component Einasto profile. 



4 Data 

Since our primary goal is to explore how well a multi-component 
DW-function describes the surface brightness profiles of ellipticals, 
we restrict ourselves to ellipticals for which a large dynamic ra- 
dial range of high resolution data is available. We hence looked at 
the well studied Virgo Cluster, for which KFCB09 provides, for 
the first time, an excellent composite compilation of many ground 
and space based observations, spanning up to five decades, i.e. the 
largest available radial range. 

In order to be able to probe the central regions, ~ 0. 1 arcsec, 
we select galaxies from the KFCB09 sample with published psf- 
deconvolved profiles for their central regions; primarily data from 
Lauer et al. (1992, 1995 & 2005). The intermediate regions of all 
galaxies are supplemented by the high resolution HST ACSVCS 
data, while the extensive outer regions come from a wide range of 
ground based observations allowing for accurate sky subtraction. 
We refer the reader to KFCB09 for details on how many datasets 
have been used for each galaxy and the details of averaging be- 
tween data sets to create the composite. 

From the sample of ellipticals in KFCB09 we exclude 
NGC4374 and NGC45I5 since they do not have psf-deconvolved 
profiles, NGC4486A since its profile is strongly contaminated by 
dust, and NGC4486B whose profile is affected by its proximity to 
M87 (NGC4486). 

NGC426I does not have psf-deconvolved profile, but we in- 
clude this galaxy in our sample. This is because a) KFCB09 uses 
HST NICMOS I.6/i H-band data for its central regions (trans- 
formed to V-magn itudes) to account for dust absorptio n in V; b) 
the 1.6/i images in lOuillen. Bower ^ Stritzingej lioool) show that 
dust absorption in near-IR is weak and restricted to the very cen- 
tral regions, and the V-H colour image shows that any nuclear point 
source is shielded by the dust; c) it has a shallow cusp and the ef- 
fects of a psf are not as strong as in steep-cusp galaxies; and d) 
during fitting with a multi-component DW-function, discussed in 
sections [STTI and l6l we varied the lower end of the fit range from 
0.07 to 0.3 arcsec (4x0.075 arcsec, the NICMOS pixel scale) and 
found that our best fit parameters were robust within 10 per cent. 

We thus obtain a sample of 23 ellipticals in and around Virgo, 
comprising of 22 ellipticals from the data set presented in KFCB0 9 
and a composite profile of NGC4494 from lNapolitano et aT] ( l2009h . 



4.1 Uncertainties in the dataset 

Critical to our modelling of multiple components is the requirement 
that the fit-residuals be consistent with measurement errors. It is 
therefore important to discuss the various uncertainties reported by 
KFCB09 for their dataset. 

Zero-points in KFCB09 are reported to have systematic un- 
certainties ^ 0.05 mag arcsec~^, and random errors of ~ 0.03 
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mag arcsec~^. While the authors do not provide quantitative val- 
ues for profile measurement errors (either for every data point or an 
rms for each galaxy) they do state that fits with a Sersic function 
are considered good when the resulting rms of fit is comparable to 
profile measurement errors, which are of the order of a "few hun- 
dredths of a mag arcsec~^ ". KFCB09 also report that the median 
rms of fits with a Sersic profile, over a restricted radial range, is 
0.04 mag arcsec~^ with a dispersion of 0.02 mag arcsec~^ . We 
thus conclude that individual galaxies in the sample can have pro- 
file measurement errors in the range of (0.02-0.06) mag arcsec"^, 
with an rms of random errors ~ 0.03 mag arcsec~ ^. 

For the central most regions ILauer et al. I l ll998f) report psf de- 
convolution errors of around 0.07 mag arcsec~^, which along 
with the random errors of 0.03 mag arcsec~^ imply that the 
central-most data points, at <0. 1 arcsec, for our selected sample 
may have uncertainties ~ 0. 1 mag arcsec^^. KFCB09 also report 
uncertainties in sky subtraction and errors due to matching profiles 
at large radii are around 0.1 mag arcsec~^. 



5 Surface Brightness with the DW-function 

In lDhar & Williamsl ( 1201 ol) the authors derive an extremely accu- 
rate analytical approximation for the line-of-sight 2D projection of 
the 3D Einasto profile, equation l ll.2t . In terms of a scaled radius 
X=-S- the surface density at a projected radial distance R is 



TiDw{R) = 



[n,b{^2X)^ 




e (5.1) 



Sb"Xe 



where. 



So = S(0) 



2eV,psnr(n) 



(5.2) 



with r(n) as the Gamma function, r[n, x] and 7(71, x] as the upper 
and lower incomplete gamma functions, respectively, and b=b{n) 
is related to rg. The parameters C2, £ , S, and ^ are functions of n, 
as derived in iDhar & WilliMnsI l l2010h (assuming b—2n and Ci = l)- 

, ,.,0 0.05657 0.00903 
C2 = 1.1513 + 5 — (5.3) 



C2 + C1 



5(X) = (C2-Ci){l"exp [-X^]} 



with, 



1.5096 0.82505 
+ ^ 



0.66299 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



As described in Paper III (in preparation), generalizing the 5 term 
( 15. 6t allows equations ( |5.U -l [5Jt to describe the surface density of 
a projected Einasto profile for any choice of scale radius r^ and an 
associated b(n). However, these equations simplify considerably 
for the N-F04 parametrization of b=2n (DWIO). This parametriza- 
tion also has another nice feature that in a log-log plot of the surface 



density, the slope of the profile is —1 at R ~ r_2 which is easily 
identified visually as the point where the profile begins to flatten out 
in log-space. In this paper we therefore adopt b = 2n and rs=r_2 
and use the parametrizations derived in Paper III to infer the 3D 
half-light radius (rsB, valid for n>0. 17) and 2D half-light radius 
(R2E, valid for 7i>0.5) of a single component through 



r-iE = r-2 



, , 0.1665 0.0035 \" 
1-5 + ^7^ 



R2E = r-2 1.5 



0.6328 _^ 0.145 



(5.8) 



(5.9) 



It should be evident from geometry that R2E<r-iE- 

Instead of fitting to estimate r-_2 and Eq, it may be more desir- 
able in some cases to fit directly for either R2E or r^E and the cor- 
responding surface brightness at Ti{R2e) and E(r3£;). In that case 
one will need to use the appropriate parametrizations for b—b{n) 
given in Paper III, which presents a number of useful results on the 
Einasto profile, and a better parametrization to measure the half- 
light radius of the 2D Sersic profile. 

The above set of equations, ( I5.1t -( |5J] |. were derived and 
tested in DWIO for 1 ^ n ^ 10. The fractional residuals with re- 
spect to a numerically projected Einasto profile were very small, 
~ 10^^ to 10^'^. In order to model the light of ellipticals, we com- 
pared equation ( 15. It to the numerical projection of il.2\ for a wider 
range of n. For 10 ^ n ^ 50 the residuals are ~ 10-^,and~ 10"^ 
for 0.5 ^ n ^ 1. The function works very well even for 0.2 ^ n 
< 0.5 with larger residuals in a small region around r^, but with an 
overall rms ^ 0.05. However, none of the galaxies require an n in 
this latter range. We also found that for the entire range 0.2 ^ n ^ 
50, the uncertainties in recovering the 3D Einasto profile parame- 
ters from 2D fits with ( 15. Il l to a numerical projection of an Einasto 
profile are better than 10^''. 

The tests were conducted over a large dynamic range in radii, 
corresponding to the domain within which the surface density drops 
from Eo to ~ 10^* Eq, which translates to a difference in magni- 
tude Afi=20 mag arcsec^^. The set of equations ( l5.U -( l577t can 
thus be used to estimate the 3D parameters of a projected Einasto 
profile for a very wide range in radii, and shape parameter 0.2 ^ n 
^ 50, corresponding to 0.02^ (a=l/n) ^ 5. 



5.1 Multi-component fits with the DW-function 

As with other 3-parameter functions discussed in section (|2j, a sin- 
gle DW-function could not fit the SB profiles over the entire dy- 
namic radial range of our galaxies. Hence, we start with the as- 
sumption that a minimum of two DW-functions, each with three 
parameters, are required to describe the light of ellipticals. 

The decision to add a third component must be based on the 
level of measurement errors in the data (section [4. It . as well as 
on the available degrees of freedom, i.e. the addition of the third 
component must be statistically justifiable. 

We decide whether to fit a 3-component model after consider- 
ing the following factors. 

(i) Overall rms of residuals: The random errors in zero-points 
of the KFCB09 data are ~±0.03ma,g arcsec~^, and result primar- 
ily from matching profiles of different filter magnitude systems. If 
the rms of residuals are much less than this level, it may indicate 
over-fitting. In that case 3-component models are not considered. 
On the other hand, if a 2-component model has a larger rms, we 
explore a 3-component model. 
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(ii) Examination of the fit residuals: A low rms of residu- 
als does not necessarily imply a good or reliable model. One can 
obtain a low rms due to very small residuals over a large radial 
range, and large residuals over a small range. We therefore exam- 
ine 2-component model residuals over the entire radial range and 
consider the fit to be good if it has consistently low residuals (< 
0.05 mag arcsec~^), except possibly at very large R, and at the 
smallest i? < 0.1 arcsec ~ 2 HST-WFPC2 pixels, where residuals 
up to 0.1 mag arcsec^^ are considered acceptable (section 143] ). 
If the overall rms is low, but there are systematically high residu- 
als in some regions, other than at very large and very small ii, we 
explore a 3-component model. 

(iii) F-test: Having fit a 3-component model, we employ the 
F-test to ensure that the reduction in rms is statistically significant 
at > 3(7, or 99.7 per cent, and is not merely due to an increase in 
the number of parameters. 

Note that a failure of the F-test does not indicate that an extra 
component is certainly not present. It just means that the number of 
degrees of freedom do not justify a statistically significant detec- 
tion of the extra component. Better resolution, observations with a 
different filter and an increase in the number of independent data 
points may lead to a significant detection of an additional compo- 
nent from the surface brightness profile. 

The resulting surface brightness is then given as: 

JV 

= ^S,(i?) (5.10) 

1=1 

with = 2 or 3 indicating the number of components and Ei (R) 
as the density of the i"* component given by equation ( 15.1b - 
i5J\ , with each component uniquely characterized by the set {Soi, 
rsi—r-2i, rii}. As indicated later in section[6jT] some galaxies may 
have A'^>3 but we do not explore this option in this paper. 

For a 2-component model we shall refer to a central compo- 
nent, which typically describes the region within the break radius, 
and an outer component, which typically identifies the main body 
of the galaxy. For a 3-component model there is an additional in- 
termediate component, which describes a transition region between 
the central and outer components, except in two cases, NGC4621 
and NGC4434, where this component indicates the presence of a 
weak system embedded within the outer component. 

It should be noted that in our models of a galaxy as a 
superposition of components, the central and intermediate DW- 
components are in excess to an inward extrapolation of the outer 
DW-component. 



6 Results of fits 

Figures rU2 1 1 show fits to the SB profiles of the 23 Virgo ellipti- 
cals with a multi-component DW-function along with the best-fit 
residuals. For all galaxies we present a 2-component model. For 
fourteen galaxies a 3-component model could be justified, while a 
2-component model is sufficient for the other nine galaxies. The 
figure captions identify whether the 2- or 3-component model is 
statistically significant. The results of fits are summarized in Ta- 
bleS 

The residuals of our models are consistent with measurement 
enws (section |4Tt over large dynamic ranges ~ 10'' in radius for 
the largest shallow-cusp galaxies down to the resolution limit of 
the HST and ~ 10® in surface brightness for the smaller steep-cusp 
galaxies. The rms is often as low as ~ 0.025 mag arcsec~^, with 
a median sample rms of 0.032 mag arcsec~^. 



The multi-component fits were carried out through a non- 
linear least squares Levenberg-Marquardt minimization using 
GNUPLOT. During the fitting process all components were allowed 
the entire dynamic radial range and no pre-defined restricted range 
in R was imposed. While convergence does depend on a reasonable 
initial guess in any non-linear fitting, we did not find any strong de- 
generacies between the fit parameters especially for models where 
residuals are consistent with measurement errors. 

As noted in section [STT] we seek a 3-component model when 
the rms of our 2-component model is greater than the 0.03 
mag arcsec"^ rms of random errors in our sample. We then per- 
form an F-test to either accept or reject the 3-component model at 
> 3(7, or 99.7 per cent. However, in the figures we present two 
cases where we do not rely on the F-test alone, and consider other 
factors. 

(i) In NGC4649 (FiglSj an F-test indicates there is a 27 per 
cent chance that the reduction in rms due to a 3-component model 
is not statistically significant; hence for the rest of the paper we 
use the 2-component model. However, from our understanding of 
errors in the central regions (section l4Tt it can be seen that the 3- 
component model certainly improves the fit near the centre and is 
acceptable at 1.2 C7. Hence, a failure of the F-test does not neces- 
sarily mean that a physically distinct system is not present. Its ex- 
istence can be verified with more information, say spectroscopic, 
about the central region. In section fTO.il we will show that the 3- 
component model may be necessary to infer the intrinsic 3D den- 
sity. 

(ii) We show a 3-component model for NGC4636 (Fig[9]( as 
an illustration where an F-test does not reject the 3-component 
model (at > 3(7), but we do. This galaxy has a large dynamic 
radial range and a 3-component model may well be admissi- 
ble. We however reject it since our first criteria to admit a 3- 
component model is that the 2-component model must have an 
rms > 0.03 mag arcsec~^ while its 2-component rms is 0.029 
mag arcsec^^. Hence, we do not feel confident in accepting this 
3-component model and emphasize that reliance on statistics and 
physical interpretations of models must be made with respect to 
the level of measurement errors in data. 

6.1 Special cases 

In this section we detail peculiarities observed while modelling 
some of the galaxies and in section 16.21 we discuss three galax- 
ies where interpretation of their components may require detailed 
modelling using additional (for example, spectroscopic) informa- 
tion. 

a) For some of the larger galaxies like NGC4486 (M87, 
Fig|2l(, even a 3-component model leaves systematic patterns ~ 0. 1 
mag arcsec~^ in the fit residuals. Hence, although the fits look 
very good and rms is better than what is typically achieved in the 
literature for such a large dynamic radial range, we believe that this 
could be an indication of an underlying fourth component. We do 
not explore four components in this paper. 

b) NGC4459 (Fig|13|l has well known embedded dust features 
which are clearly evident in the SB profiles. Its outer n of the 2- 
component model is fairly robust with respect to whether we in- 
clude or exclude the radial range affected by the dust. This is, how- 
ever, not true of its inner n, which changes appreciably based on 
the inclusion or exclusion of the dusty region. We adopt the fit that 
includes the dust region. 

c) NGC4473 (Fig|14b is quite an interesting case. The SB is 
fit extremely well with an rms of 0.031 mag arcsec^^ using only 
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Table 1. Multi -component DW model properties of Virgo Ellipticals 



Component Parameters Effective Radii RMS (mag) 

Name Type D Scale Uc ntnt 'lo fuEc f'zEint f'zEo f'zE R2E /^o Qc Av Vt Myt Lyt Multi- Double- 

(Mpc) (pc) (pc) (kpc) (kpc) (kpc) (kpc) (mag^) (mag) (mag) (xIO^Lq) DW Sersic 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (II) (12) (13) (14) (15) (16) (17) (18) (19) (20) 



NGC4472 


E2 


17.14 


83.10 


0.717 


1.860 


5.359 


382.75 


2.09 


29.87 


25.59 


18.88 


15.92 


0.806 


0.072 


7.88 


-23.29 


176.5 


0.031 


0.07 


NGC4486 


El 


17.22 


83.48 


1.090 


2.619 


6.451 


30.81 


3.96 


75.26 


50.98 


37.29 


16.21 


0.722 


0.072 


8.04 


-23.14 


153.9 


0.050 


0.09 


NGC4649 


E2 


17.30 


83.87 


0.997 




5.444 


951.50 




14.96 


13.78 


10.20 


15.64 


0.828 


0.086 


8.46 


-22.73 


106.2 


0.032 


0.08 


NOC4406 


E3 


16.83 


81.59 


0.895 


1.488 


3.084 


240.85 


32.58 


4.14 


24.61 


17.96 


15.66 


0.709 


0.096 


8.39 


-22.74 


106.9 


0.049 


0.13 


NGC4365 


E3 


23.33 


113.11 


0.970 


3.720 


7.233 


549.05 


5.07 


43.98 


28.07 


20.70 


15.58 


0.717 


0.068 


9.09 


-22.75 


107.5 


0.036 


0.09 


NGC4261 


E2 


31.60 


153.20 


0.896 


1.446 


5.023 


533.37 


2.35 


22.45 


16.29 


11.92 


16.11 


0.794 


0.059 


9.90 


-22.60 


93.45 


0.024 


0.08 


NGC4382 


E2 


17.86 


86.59 


1.545 


0.903 


3.330 


186.71 


0.70 


12.21 


11.14 


8.25 


14.65 


0.761 


0.101 


8.78 


-22.48 


83.74 


0.090 


0.11 


NGC4636 


E3 


14.70 


71.27 


1.282 




6.064 


632.50 




26.25 


25.44 


18.85 


16.41 


0.760 


0.090 


8.54 


-22.30 


70.82 


0.029 


0.04 


NGC4552 


El 


15.85 


76.84 


0.754 


3.306 


6.528 


89.56 


1.09 


20.34 


12.51 


9.14 


14.48 


0.873 


0.133 


9.29 


-21.71 


41.22 


0.049 


0.09 


NGC4621 


E4 


14.93 


72.38 


3.714 


0.985 


9.561 


16.13 


4.67 


8.77 


7.62 


5.72 


10.21 


0.742 


0.107 


9.29 


-21.58 


36.75 


0.027 


0.05 


tNOC4494 


El 


15.85 


76.84 


1.701 


2.057 


4.171 


10.95 


0.19 


5.54 


5.15 


3.81 


12.61 


0.838 


0.067 


9.90 


-21.10 


23.50 


0.025 




NGC4459 


E2 


16.07 


77.91 


4.091 




3.835 


257.37 




4.84 


4.44 


3.29 


13.16 


0.804 


0.149 


10.09 


-20.94 


20.36 


0.050 


0.05 


NGC4473 


E4 


15.28 


74.08 


2.300 




5.649 


532.50 




5.48 


4.26 


3.13 


14.44 


0.607 


0.092 


10.00 


-20.92 


20.03 


0.031 


0.05 


NGC4478 


E2 


16.98 


82.32 


0.941 


1.217 


2.641 


4.25 


0.07 


1.49 


1.45 


1.08 


13.55 


0.822 


0.080 


11.37 


-19.78 


6.99 


0.035 


0.10 


NGC4434 


EO 


22.39 


108.55 


1.238 


0.576 


5.532 


15.50 


2.78 


1.47 


1.73 


1.27 


13.61 


0.928 


0.072 


12.18 


-19.57 


5.74 


0.032 


0.07 


NGC4387 


E4 


17.95 


87.02 


3.562 




2.621 


82.31 




1.68 


1.65 


1.23 


14.18 


0.633 


0.107 


12.14 


-19.13 


3.83 


0.037 




NGC4551 


E3 


16.14 


78.25 


1.634 


1.214 


2.446 


23.12 


0.15 


1.67 


1.59 


1.18 


14.47 


0.734 


0.125 


11.95 


-19.09 


3.69 


0.026 


0.05 


NGC4458 


El 


16.37 


79.36 


1.502 


3.078 


3.149 


23.05 


0.12 


2.07 


1.74 


1.28 


12.93 


0.879 


0.077 


12.12 


-18.95 


3.24 


0.026 


0.06 


NGC4464 


E3 


15.85 


76.84 


1.994 


1.222 


3.094 


13.66 


0.08 


0.78 


0.69 


0.51 


12.74 


0.749 


0.071 


12.59 


-18.41 


1.97 


0.021 


0.06 


NGC4467 


E3 


16.53 


80.14 


2.673 




2.402 


31.67 




0.54 


0.51 


0.38 


14.38 


0.813 


0.074 


14.18 


-16.92 


0.499 


0.020 


0.02 


VCC1440 


EO 


16.00 


77.57 


1.378 




4.976 


8.79 




0.91 


0.90 


0.67 


14.36 


0.965 


0.088 


14.14 


-16.88 


0.484 


0.032 


0.04 


VCC1627 


EO 


15.63 


75.78 


2.002 




2.907 


17.47 




0.38 


0.37 


0.27 


14.63 


0.928 


0.127 


14.53 


-16.44 


0.324 


0.037 


0.04 


VCC1199 


El 


16.53 


80.14 


1.896 




2.389 


13.49 




0.23 


0.21 


0.16 


14.22 


0.869 


0.071 


15.55 


-15.54 


0.140 


0.023 


0.02 



NOTES. — Galaxy prop erties from the multi-component DW models best-fitting the surface brightness (SB) profiles coinpiled from KFCB09, except that of NGC4494 (marked with a |) whose data is from 
iNaoolitano et alj {200^: we have, however, transformed their intermediate-axis profile to major-axis so as to have uniformity with the KFCB09 sample. Columns contain: (2) Galaxy type as defined in KFCB09. 
(3) Distance in Mpc from KFCB09 and references therein. (4) shows the physical size of 1 arcsec in pc at the distance given in column 3. (5-10) Einasto shape pai'ameter n and the intrinsic (3D) effective or 
half-light radius r^j^ of the central-, intermediate- (when applicable) and outer components, deduced from the multi-component DW fits to the SB profile, r^^j^ of the components have been computed from the 
respective best fitting r-2 using equation 15. 8L (1 1-12) Total effective or half-light radii of the galaxy, intrinsic (vse) and projected (i?.2£;). deduced by numerically integrating the best fitting multi-component DW 
inodel to infinity. (13) V-band Central (r— 0) suiface brightness in mag arcsec~^ deduced from the best fitting multi-component DW model and corrected for Galactic extinction (column 15). (14) characteristic 
axis-ratio of the galaxy (see section l6Ji . (16-17) Galactic extinction corrected total V-band apparent (Vr) and absolute {Mvt) magnitude. (18) V-band extinction corrected luminosity assuming A/^q— 4.83, 
(19-20) Comparison of rms of residuals between the best-fitting multi-component DW models (this paper) and double-Sersic models in Hopkins et al. (2009a,b), for the same SB dataset of KFCB09. 



a 2-component DW-function. The outer component has an Einasto 
shape parameter n=5.65 and the central region shows a density 
profile flattening typical of shallow-cusp galaxies. It has A'Ivt=- 
20.92, Lv-T=2xlO^"i0 and rather elliptical (E4/E5) isophotes, 
typical of the steep-cusp galaxies. The KFCB09 composite data 
extends to 300 arcsec (22 kpc), an extent that is larger than that 
of all steep-cusp ellipticals (except the unusual NGC4621, section 
I6.2t . but smaller than that of the smallest shallow-cusp elliptical 
NGC4552 (FigllOt. It has a Galactic extinction corrected central 
surface brightness of 14.44 mag arcsec~^ which is also found in 
both families of galaxies (Tabl^ However, its intrinsic 3D central 
density as well as the overall intrinsic density profile appears to be 
similar to that of shallow-cusp galaxies. 

d) NGC4387 (Fig[T4j is well fit with a 2-component model but 
may well have a third component. Since the 2-component rms is 
already reasonably low at 0.0367, we do not explore a 3-component 
model. However the SB profile is very similar to that of NGC4551 
(FiglITt which clearly shows a 3-component structure. The central 
component for NGC4387 is also spatially more extended than that 
of NGC455 1 and has a larger n compared to its outer component; a 
feature that is different from most galaxies in the sample. The n and 
r_2 of the outer component for NGC4387 and NGC4551 are also 
very similar. This suggests that the central component of NGC4387 
could well be a sum of two components, which a SB analysis with 
the present data is not able to distinguish. 

Note that the fading light of the giant elliptical NGC4406 
(M86) beyond 575 arcsec (Fig|4]and section |6^ . affects the en- 
tire region of NGC4387, The centres of these two galaxies are sep- 
arated by 668 arcsec. Their SB profiles extend through 800 and 
93 arcsec for NGC4406 and NGC4387, respectively. The high SB 



central regions of NGC4387 are unlikely to be affected by the low 
signal (> 25 mag arcsec^^) from NGC4406, but the outer regions 
beyond ~ 50 arcsec could be. However, the SB of NGC4387 show 
no detectable features in its outer profile, 

e) VCC1440 (Fig|20b is unusual in that despite being a fairly 
low luminosity steep-cusp galaxy (Mvt = —16,88, I/v't=3,2x 
10* Lq) it has a fairly large 7i=4,98 for its outer component, sim- 
ilar to the massive shallow-cusp galaxies. 



6.2 Galaxies with uncertain components 

The analysis of the SB profiles of NGC4406, NGC4382 and 
NGC4621 do not lead us to an unambiguous conclusion of which 
component characterizes the main body of the galaxy. We hence ex- 
clude these galaxies from studies involving specific aspects of indi- 
vidual components, for example, trends involving shape parameters 
and luminosity of components. The figure captions identify them as 
exceptions. The details are explained below: 

a) The SB of NGC4406 (FigEJ is fit well with a 3-component 
model with an rms of 0,049 mag arcsec^^. However, it is not clear 
what constitutes the main body of the galaxy. The 7i=3,084 compo- 
nent contains 23 per cent of the total light, while the outer compo- 
nent contains 76 per cent, but has an unusually low n= 1.488. Note 
that while NGC4406 has overlapping profiles with NGC4387 and 
NGC4374, they overlap only beyond 550 arcsec. The third com- 
ponent gains prominence around 100 arcsec and continues through 
800 arcsec. It is hence intrinsic to NGC4406 and not a feature due 
to incorrect subtraction of the light of NGC4387 and NGC4374, 

This galaxy is streaming into Virgo at 1400 km s^'^ and it 
is possible that it has also gone through a recent interaction or a 
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merger. Chandra images and the presence of large plumes of HI 
gas around this galaxy do indicate such a possibility. Therefore it 
is not surprising that the n of its outer component is different from 
the typical values of other galaxies. Same is true of the intermediate 
component. 

Our 2-component fit over four decades in radius (with an up- 
per limit of 153 arcsec, as in KFCB09) is not a good fit, although its 
outer 71 is consistent with that of other large shallow-cusp galaxies. 
However, as mentioned earlier estimating fit parameters and struc- 
tural properties from a limited radial range can be misleading. For 
example, the half-light estimated with the 2-component model is 
31 kpc, much larger than the 3-component half-light of 18 kpc. 

In Fig|4]and Figs. | |22I231 and l24t we identify the component 
with n=3.084 as the one characterizing the main body of the galaxy 
(which for most galaxies, is the outer component) and the com- 
ponent with n= 1.488 as a perturbing or transitional component 
(which for most galaxies, is the intermediate component). 

b) In the 3-component model for NGC4382 (Fig|7j, the outer 
component contains 94 per cent of the total light. However, the 
rather abnormal bump in the profile around 100 arcsec and the large 
residual patterns may indicate that the system has not relaxed since 
a recent interaction or merger event (refer to KFCB09 for an image 
from |http://www.wikisky.org| showing signatures of a recent inter- 
action). 

Hence the outer n of its 3-component model may not be repre- 
sentative of well relaxed systems, although its central and interme- 
diate components are similar to those of other shallow-cusp galax- 
ies analyzed in this paper. In fact, if the bump around 100 arcsec 
is excluded (Fig[8]l, we find an outer n of a 2-component model 
more consistent with those of other shallow-cusp galaxies. How- 
ever, we do not want to draw inferences on components obtained 
after excluding some radial ranges from the fit, especially in the ab- 
sence of a theoretical motivation for doing so. We hence adopt the 
3-component model in Fig|7] 

c) NGC4621 (Figinil is an interesting case. With Myr^- 
21.58, it is the most luminous steep-cusp galaxy in our sample. In 
general, galaxies with Mvt<-21.5 appear to have shallow-cusps, 
and are spatially more extended, while those with Mvt>-21.5 
have much smaller spatial extent and usually have steep cusps. 
However, NGC462 1 has the steepest central cusp and has the high- 
est central surface brightness of all galaxies in our sample and is 
physically almost as large, and as luminous as the smallest shallow- 
cusp galaxy NGC4552 (Fig[lOll. 

The residuals of a 2-component fit and an F-test justifies a 3- 
component model, which also gives a far more reasonable 7i=3.71 
for its central component than the unusually large n=8.79 of a 2- 
component model. However, the outer n of its 3-component model 
is much larger than that of any galaxy in our sample, including the 
shallow-cusp galaxies that generally have large outer n. This could 
be due to the embedded intermediate component with n=0.985. 
NGC4621 is a E4 galaxy and its elliptical isophotes show sharp 
pointed features. It could be that the ellipticity of the embedded 
intermediate component is quite different from that of the outer 
component and accounting for this ellipticity may be necessary to 
obtain a better estimate of n values. At the same time this galaxy 
could just be a special case, as mentioned above. 

NGC4434 (Fig|16|( may be a similar system. Its intermediate 
component, which is nearly gaussian, with n=0.576, is embedded 
within its outer component. The galaxy however is fairly round 
(EO) with regular isophotes. We thus do not exclude this galaxy 
but note that like NGC462I and NGC4473 and VCC1440 (section 
I6.lt its outer n is larger than that of all other steep-cusp galaxies. 



and is are more consistent with that of the shallow-cusp galaxies. 
However, like NGC4621, the embedded component may have al- 
tered the shape of the outer component. Its outer n may therefore 
have large uncertainties. 

6.3 Estimating luminosity and half-light radii 

In order to calculate the total luminosity, the cumulative luminosity 
within R and the half-light radii one has to, in principle, account for 
the varying ellipticity e—1 — q; where q=b/a is the axis ratio. Nu- 
merically, this can be done by expressing the area element dA{R) 
in terms of an axis ratio, q{R), as 



dA{R) = 2Tvq{R)R + TvR' 
and the projected luminosity is then given by 



dR 



Lp{Ra) — 



E{R)dA{R) 



(6.1) 



(6.2) 



where R and Ra are along the major axis. However, this involves 
taking derivatives of axis ratios which show large and some times 
abrupt variations. Further, in addition to real variation in q there 
could also exist artificial variations due to the limitations of the 
ellipse-fitting process. 

To avoid the above difficulties, during the ellipse-fitting pro- 
cess one can add up the projected light non-parametrically in el- 
liptical isophotes, L^^ {Ra), through the last data point, Ra, and 
define a characteristic axis ratio, qc, such that 



Lp^ {Ra) 



27r/"" T,AN{R)RdR 



(6.3) 



where 'Ean{R) is the spherically symmetric analytical function; 
here a multi-component DW-function. 

To estimate the total luminosity, we use this value to integrate 
through infinity, even though this ellipticity seldom represents the 
axis ratio of the outer most isophotes. Doing so can be justified for 
profiles with large dynamic radial range, as in this paper, since by 
definition (equation i6.3i ) it is a weighted axis ratio and the ex- 
tremely low luminosity regions extending through infinity are un- 
likely to change this weight. Further, since we do not know the true 
axis ratio beyond the last data point we would not like to assume 
that the axis ratio through oo is equal to that of the last data point. 

The total luminosity is then given by 



LpT — 2-Kqc 



{R)RdR 



(6.4) 



and the projected half-light radius, R2E, can be estimated from 

/R2E 
T,AN{R)RdR. (6.5) 



Note that the characteristic axis ratio defined above, i6.3h is 
neither the mean nor a luminosity weighted axis ratio in the usual 
sense. Nevertheless, it is a useful definition that reduces the error in 
estimating the total light and half-light radii, and also ensures that 
they are weakly dependent on the specific choice of parametrization 
of the SB profile and assumptions of ellipticity; provided of course 
that the fit parameters have been deduced after modelling the entire 
dynamic radial range and not from a limited range. 

Analytical expressions for the integrals in equation ( I6.3l l and 
i6.4\ may not exist. However, for the projection of an Einasto pro- 
file, analytical forms are given in Paper III, which accurately mod- 
els a numerical integration with the DW-function as well. 



10 



Barun Kumar Dhar and Liliya L.R. Williams 




Figure 1. Multi-component DW models of the surface brightness (SB) profiles of 23 galaxies in the Virgo Cluster. For each galaxy, the figure keys 
show the best-fitting Einasto index n and scale radius rs=r_2 of the components. The rms of residual is shown in the adjoining residual profile. 
The fit range is listed in the figure panels and in a few cases when some points are excluded from the fit, the excluded region is also marked on 
the SB profile. A 2-component model is always shown, while a 3-component model is shown when it is found to be statistically significant through 
an F-test (section IsTTl . The model accepted is marked 'adopted' in the caption. Circles indicate data from KFCB09, solid (red) line shows the total 
multi-component DW model, dot-dashed (blue) marks the central component, short-dashed (black) the outer component or the component resembling 
the main body of the galaxy and long-dashed (magenta) shows an intermediate or embedded component in some galaxies. Note that the central and 
intermediate DW-components are in excess to an inward extrapolation of the outer DW-component in those regions. 

Above NGC4472 (VCC1226) Left: 2-component, Right: 3-component (adopted). The large residuals around 1.5 arcsec in the 2-component model 
indicates that a 3-component model may be necessary, which not only has a much lower rms of 0.0309 but also has consistent low and non-divergent 
residuals < 0.05 over a large dynamic radial range ~ 10^. Also note the occurrence of an Ernasto index of n < 1 in the central component of the 
statistically significant (adopted) model; as in all massive shallow-cusp galaxies. (Colour versions of these figures are available in the online edition.) 



In this paper, we use the apparent magnitudes estimated non- 
parametrically by KCFB09 through the last data point to calculate 
the numerator in equation J6.3t . As a cross-check one can verify 
that the characteristic axis ratio estimated through equation i6.3\ 
is usually consistent with the E-type of the galaxies listed in Ta- 
ble [T] and the ellipticity values in the literature. The V-band total 
(integrated through infinity) magnitudes are listed as Vt and the 
corresponding total absolute magnitudes are listed under Mvt- 
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Figure 2. NGC4486 (VCC1316 or M87) Left: 2-component; Right: 3-component (adopted). The 3-component model shows an overall improvement and a 
reasonable rms=0.0497 with non-divergent residuals over an extremely large radial range ~ 10^. However, the presence of (~ 0.1 mag) systematic patterns 
between 3-150 arcsec, greater than the ~ 0.03 mag measurement errors, suggests the possible presence of a fourth component. This is not surprising given the 
large spatial extent of this galaxy. We do not explore 4-component models in this paper. Refer to caption of Fig[T]for further details. 




Figure 3. NGC4649 (VCC1978) Left: 2-component (adopted); Right: 3-component (significant at \.2a). The 3-component model could be used to compare 
with non-parametric deprojections (sectioi fTO. U . Refer to caption of Fig[T]for details. 
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Figure 4. NGC4406 (VCC881) Left: 2-component fit upto 153.3 arcsec - the upper-limit of fit with a Sersic profile shown in KFCB09; Right: 3-component 
(adopted) fit over the entire radial range through 800 arcsec. However, which component forms the main body of the galaxy is not clear and hence this galaxy 
is an exception (section l6!2t . Refer to caption of Fig[T]and section l6!2l 




Figure 5. NGC4365 (VCC731) Left: 2-component; Right: 3-component (adopted). Refer to caption of Fig[T]for details. 
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Figure 6. NGC4261 (VCC345) Left: 2-component; Right: 3-component (adopted). The improved fit with a 3-component model, around the transition radius, 
and the near vanishing of large scale systematic patterns along with a significant reduction in rms is cleai'ly visible. Refer to caption of Fig[T]for details. 
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Figure 7. NGC4382 (VCC798) Left: 2-component; Right: 3-component (adopted). Although this galaxy shows a visibly shallow-cusp, note the lower outer-n 
and an inner n = 1.55 similar to that of some of the less massive steep-cusp galaxies. The 3-component fit is very good, except in a region around 100 arcsec. 
Hence the best-fitting parameters of its outer component may not represent the true values. The outer regions of this galaxy also show indications that it may 
not have relaxed from a recent interaction and we therefore mark this galaxy as an exception (see section |6^ . Refer to caption of Fig[T]for further details. 
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Figure 8. NGC4382 (VCC798) Left: 2-component; Right: 3-component. As noted in Fig. [T] tliis galaxy is an exception. Here we show a fit by excluding a 
region around the bump from 26 — 221 arcsec, as in KFCB09. Note the larger outer-n and the n ~ 1 for the inner component, similar to the pattern seen in 
other massive shallow-cusp galaxies. However, we adopt the 3-component fit in Fig.|2]since excluding domains of fit with ad-hoc fitting functions, may not 
reveal the true structure of these regions. Refer to caption of Fig[T]for further details. 
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Figure 9. NGC4636 (VCC1939) Left: 2-component (adopted); Right: 3-component. The 3-component model is shown as an illustration of a case where an 
F-test does not reject it but we do, since the 2-component model has a rms~0.03 and the residuals are non-divergent. Refer to caption of Fig[T]for details. 
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Figure 10. NGC4552 (VCC1632) Left: 2-component; Right: 3-component (adopted). At Mvr=— 21.71, the smallest shallow-cusp elliptical in Virgo is more 
luminous than all other steep-cusp ellipti cals. As in most shall o w-cusp ellipticals, it als o has an n<l for its central component which is not seen in any of the 
steep-cusp elhpticals. HST observations iRenzini et alj fl993) . ICaDDellari et al] il999l) ^ reveal variable UV-flare activity in the centre which is inteipreted to 
arise from a low-level AGN. Contributions from such a point source to the central-most data point is excluded from the fit. See caption of Fig[T]for details. 




Figure 11. NGC4494. Left: 2-component; Right: 3-component(adopted). This ga laxy is in the Coma-I c loud around Virgo and comparable to the more 
luminous steep-cusp (> 2 X IO^^L^q) Virgo ellipticals. The composite data from lNapolitano et al.1 )200^ is plotted in terms of the major axis radius, as in 
other galaxies from the KFCB09 sample. This is the only galaxy that is not in the KFCB09 sample. Refer to caption of Fig[T]for details. 
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Figure 12. NGC4621 (VCC1903) Left: 2-component; Right: 3-component (adopted). The most luminous (Mvt=-21.58, Lvt=3.68x lO^^Lv-©) steep- 
cusp elhptical in Virgo. The 3-component model improves the fit in the region from 9-130 arcsec and also gives a more reasonable central-n. However, it 
also has a much larger outer-n. The SB has shaip pointed isophotes and it is not clear how much the embedded component and its ellipticity influences our 
determination of outer-n. We hence mark this galaxy as an exception (see section l6!2l . Also see caption of Fig[T] 




Figure 13. NGC4459 (VCC1154) Left: 2-component(adopted); Right: 2-component excluding the region affected by dust from 1.24 — 9.6 arcsec, as in 
KFCB09. It is one of three galaxies in our sample that has an rJcenira! > "-outer. Refer to caption of Fig[T]for details. 
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Figure 14. Left NGC4473 (VCC1231): Its 2-component model has an mis of only 0.031 over a dynamic radial range ~10'' and 12 mags in SB. It is the least 
luminous of r^W^'^ LyQ ellipticals in Virgo with properties similar to both the steep and shallow cusp families (section ISTt . However, unlike the massive 
shallow cusps that have central-n<l, it has n=2.23 leading to a much gradual central flattening than the sharper transition seen in the former (section [Ot ; 
Right NGC4487 (VCC828): It is one of three steep-cusp ellipticals with ni.entral>i^'oufer and has a luminosity of 3.83x IO^L^^q. Its profile could be 
affected by the light of NGC4406 (section lSTl . Note the similarity of this profile with that of NGC4551 (Fig. II 71 where 3-components could be justified. 




Figure 15. NGC4478 (VCC1279) Left: 2-component; Right: 3-component (adopted). As in NGC4494, three components are clearly visible. It is the most 
luminous of ~10^L^/q ellipticals in Virgo. Along with NGC4434, they are the only two steep-cusp ellipticals with a central-n~l. See Fig[T]for details. 



1 8 Barun Kumar Dhar and Liliya L.R. Williams 




Figure 16. NGC4434 (VCC1025) Left: 2-component; Right: 3-component (adopted). A 3-component model is clearly needed. However, it has an unusually 
large outer-n typical of the more massive shallow-cusp galaxies. Refer to caption of Fig[T]for details. 




Figure 17. NGC4551 (VCC1630) Left: 2-component; Right: 3-component (adopted). Refer to caption of Fig[T]for details. 
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Figure 18. NGC4458 (VCCl 146) Left: 2-component; Right: 3-component (adopted). Refer to caption of Fig[T]for details. 
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Figure 19. NGC4464 (VCCl 178) Left: 2-component; Right: 3-component (adopted). It is the least luminous of ~ lO^LyQ ellipticals in Virgo. Refer to 
caption of Fig[T]for details. 
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Figure 20. Left NGC4467 (VCC1192); Right: VCC1440. Low luminosity ellipticals in Virgo (~ lO**Ly0) bordering the dwarf elliptical population. Both 
galaxies show a steep-cusp, but note the large outer-n for VCC1440 typical of the more massive shallow cusps. Refer to caption of Fig[T]for details. 
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Figure 21. Left: VCC1627; Right: VCC1199. The least luminous elHpticals in Virgo (~ 10«L^q). Refer to caption of Fig[T]for details. 
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7 Comparison with other parametrizations 

The goal of this paper is to investigate how well can a multi- 
component DW-function model the SB and consequently whether 
the intrinsic 3D density can be described with a multi-component 
Einasto model. A detailed comparison with other parametrizations 
for every galaxy is beyond the scope of this paper. However, these 
parametrizations often show residuals larger than measurement er- 
rors (sections l3.1l and [T2l > . Since the literature has a few examples 
of fits with other parametrizations, for the galaxies we model here, 
we present a comparison in sections l7!T]|7.2l and l7.3l 



7.1 Core-Sersic, SersicH-King and Nuker models 

All galaxies in our sample have been fit with a combination of 
psf-convolved Sersic, Core-Sersic (hereafter, CS) and King mod- 
els to the ACSVCS Sloan g- and z-band profiles in lperrarese et al.l 
( l2006l) . While the residuals are difficult to ascertain from their fig- 
ures, fit residuals of King-l-Sersic, Core-Sersic and Nuker models in 
th e central 0.02 — 20 arcsec, or ~ 3 radial decades, are presented 
in lLauer et al] j2007h (L-l-07); see their fig. 9 and 10. These are the 
same psf-deconvolved profiles that were used in KFCB09 and the 
present paper where the data in the central regions are from psf- 
deconvolved WFPC V-band and 1.6 jj, NICMOS images of Lauer 
et.al. (1992, 1995 & 2005). Note that L-l-07 fit the same models used 
in F-l-06 (for a given galaxy) to psf-deconvolved WFPC/NICMOS 
images and then compare the F-l-06 models to fits with a Nuker 
model. L-l-07 also notes that differences in choice of psf, camera 
and observing band are not significant to prevent a comparison. 

The Nuker fits usually exclude fitting 'nuclei'. These regions 
are fit by F-l-06 with a King and a Sersic (6-parameters) or CS pro- 
file (9-parameters); the CS-l-King galaxies of F-l-06 are not shown 
in L-l-07. Further, the fits in L-l-07 are compared within the cen- 
tral 10 — 20 arcsec domain of validity of the Nuker profile. All 
other parametrizations have been fit over larger radial ranges. Fits 
in this paper, which use KFCB09 data, use the largest range. The 
19 galaxies that overlap with our sample are discussed below: 

(i) Steep-cusp with a King model for 'nuclei' : For most steep- 
cusp galaxies F-l-06 required a King profile for the central region 
and identified them as 'nuclei'. Note that King models have a 
truncation radius and a flat core in 3D, and hence in 2D as well. 
The Nuker fits exclude these regions. For galaxies that we have in 
common with F-l-06 and L-^07 - NGC4387, NGC4467, NGC4551, 
NGC4458, VCC1199, VCC1440 and VCC1627 - we show that 
they can be well fit with DW-functions, often with just two com- 
ponents over the entire dynamic radial range, and with much better 
residuals than with either the King-l-Sersic or the Nuker profiles. 

(ii) Steep-cusp with a single Sersic model: For some galaxies, 
F-l-06 showed that the entire profile can be described with a sin- 
gle component Sersic model. The comparison plots in L-l-07 show 
that a single Sersic profile can not fit the central-most regions of 
NGC4478, NGC4473, NGC4621, NGC4434 and NGC4464 where 
the Nuker performs better. We show that it is possible to easily 
quantify these regions using the DW-function with better residuals 
than with the Nuker profile. 

(iii) Shallow-cusp: NGC4365, NGC4382, NGC4406, 
NGC4552, NGC4649, NGC4472 and NGC4486 have been fit 
with the Core-Sersic and Nuker profiles. Although the residuals 
are smaller than those in the case of the steep-cusp galaxies, we 
show that fits with a multi-component DW-function produce even 
smaller residuals. It should also be noted that the Nuker fits apply 



to a limited range, and as shown by iGraham et alj l l2003h the 
parameters depend strongly on the selected radial domain of fit. 

7.2 Single-Sersic models 

KFCB09, whose composite SB profiles we use in this paper, 
presents fits to the SB with a single Sersic 3-parameter function for 
each galaxy. Since such a 1 -component Sersic profile cannot ade- 
quately fit the full radial range of 4-5 decades of KFCB09 profiles, 
the authors estimate the largest radial range over which a single 
Sersic profile produces robust fits and residuals comparable to mea- 
surement errors. This range is typically 1-2.5 radial decades, and 
excludes regions interior to the transition radius, where the slope 
changes rather abruptly. 

KFCB09 also had to invoke certain constraints to estimate the 
total luminosity and half-light radii, for example, a limiting magni- 
tude up to which to integrate the light. This was usually done for 
galaxies where the Sersic profile failed to model the SB profile over 
a large radial range. While as noted in KFCB09 the existence of a 
physically justified limiting magnitude is possible - for example, 
a tidal truncation radius or incorrectly subtracted light of a neigh- 
bour - such limits should not exist due to the failure of an ad-hoc 
fitting function to model the entire SB profile. Further, since the fits 
were obtained from a limited range, they give a biased estimate of 
the Sersic index (for example, in the case of NGC4406 and other 
massive galaxies in table 1 of KFCB09). 

We agree with KFCB09 that fits must be consistent with the 
measurement errors, and show that this can be achieved using a 
multi-component DW-function, which also avoids systematic de- 
viations between the data and the fits. Our estimates of structural 
properties are therefore likely to be more meaningful, and can be 
used to obtain direct estimates of the intrinsic 3D structural proper- 
ties of galaxies. 

7.3 Double-Sersic models 

Hopkins et al. (2009a,b) fit a 2-component Sersic profile to all 
galaxies in the KFCB09 sample that we use here. While they do 
not show the residual profiles, we note that the rms of their fits 
with a 2-component Sersic profile is usually 40 per cent larger than 
the rms of our 2-component DW models, and nearly 200 — 300 
per cent larger than the rms of our 3-component DW models. In 
Table [T] we list the rms of fits of our best-fitting DW models and 
that of the double-Sersic models in Hopkins et al., for comparison. 

From their double Sersic models Hopkins et.al. conclude that 
steep-cusp and shallow-cusp galaxies are two disjoint populations, 
and that the outer Sersic index m does not depend on the mass or 
luminosity of these galaxies. 

Our fits lead us to a different conclusion. We observe that the 
outer Einasto index n does increase with luminosity (and conse- 
quently size, rsE - the 3D half-light radius) in a s eemingly con- 
tinuou s manner. A similar trend was also noted by iGraham et al.l 
( 1 19961) and F-l-06 based on their fits with Sersic profiles. 

7.4 The form of the 2D structure 

An examination of the SB profiles show that the central regions 
have distinct variations in slope. This is more true for the steep-cusp 
galaxies than shallow-cusp ones. This deviation from pure power- 
laws are reflected in the large fit residuals in the central regions as 
discussed in section FTTT] 
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Fig. lU2ll when compared to the fits in the literature described 
above, demonstrate that the 2- or 3-component DW-function pro- 
vides a better fit, over a larger radial range, than other existing func- 
tional forms. Our overall r?ns are comparable to, or lower than that 
of other models, and our residuals are consistently low over the 4-5 
radial decades of the available composite observations (KFCB09). 
Furthermore, the DW-function is a very accurate 2D projection of 
the 3D Einasto profile, and is expressed in terms of the 3D Einasto 
profile parameters. This means that if the 2D fits are good, the in- 
trinsic 3D luminosity structure is that of superimposed Einasto pro- 
files, and can be inferred directly, with no further modelling. 

We hence propose that the light of ellipticals that was believed 
to be well fit with a Sersic profile in 2D, is instead better described 
by a multi-component form of a similar function (the Einasto pro- 
file) in 3D, whose 2D projection is given by a multi-component 
DW-function. 



8 Component properties: Luminosity, half- 
light radius and Einasto index n 

In this section we investigate the structural properties of the compo- 
nents deduced from the multi-component DW fits. We shall be re- 
ferring to the statistically significant best-fitting models only, listed 
in Table[T] and not all 2-component and 3-component fits shown in 
Figs lll2l1 Nine galaxies are described with two DW-components 
(two shallow cusps and seven steep cusps), and fourteen galaxies 
are described with three DW-components (seven each of shallow 
and steep cusp galaxies). 

In our modelling of a galaxy as a linear superposition of com- 
ponents, the central and intermediate DW-components described in 
the following sections 18.118.31 are in excess to an iimer extrapola- 
tion of the outer DW-component; and do not contain all of the light 
in the central and intermediate regions. This is consistent with sim- 
ilar decompositions in the literature. Superposition of components, 
comprising of a Sersic profile for the central bulge superimposed 
on an underlying exponential model (a Sersic profile wi th m = 1), 
is ofte n applied to the case of lenticulars and spirals. ICote et al] 
( |2007|) have shown that a similar decomposition using a double- 
Sersic profile as a fitting function can also be applied to ellipticals 
with Mb > —19.5. The resulting central Sersic-component is then 
used to evaluate physical properties of the central region. 

For the case of shallow-cusp gala xies, as in the spectro- 
scopic and kinematic modelling of M87 in lTenies. Einasto & HaudI 
( Il99lh and the fitt ing function (double S ersic) based modelling of 
'core' galaxies in iHopkins et al] l l2009bh . we also find that these 
galaxies can be modelled as a linear superposition of components. 
While such fit components need not correspond to real physi- 
cal systems, we show in section 18.41 three cases of shallow-cusp 
('core') galaxies whose central DW-components coincide very well 
with spectroscopically identified systems. 

In order to estimate luminosities of components, a knowledge 
of a characteristic axis ratio (equation i6.3\ ) over the extent of each 
component is required. Since the components are superimposed, it 
is very difficult to isolate their characteristic axis ratios using SB 
analysis; except maybe the outer-most dominant one. However, el- 
liptical galaxies have axis ratios of 0.3 (E7) ^ g ^ 1.0 (EO) and for 
any arbitrary choice of q the uncertainty in estimating the luminos- 
ity can at most be off by a factor of 2 — 3. In our sample, except 3 
galaxies which are E4, the rest are between El and E3. Our worst 
errors in estimating component luminosities are therefore less than 



a factor of 2. Additional spectroscopic or kinematic modelling may 
be used to constrain the ellipticities of individual components. 

With this understanding, for the purpose of computing lumi- 
nosity and half-light radius of components, we assume that the axis 
ratio of all components are the same as the characteristic axis ratio 
Qc for the entire galaxy (equation J6.3b ). Note that uncertainties in 
interpreting component luminosities do not affect our estimate of 
the total luminosity of the galaxy. 



8.1 Central component 

From the nine galaxies with two components and fourteen galaxies 
with three components, we observe that: 

1) In all galaxies, the central component has a lower n than 
that of the outer component. Since Einasto functions of the form 

_ 2 _ 

e ^ (n=0.5) — >■ faster than e ^ (n—l), a lower n implies a 
larger concentration. Hence, the central component is more con- 
centrated than the outer component. 

2) The central component of the shallow-cusp galaxies are 
one or two orders of magnitude more luminous than that of the 
steep-cusp galaxies (Fig|22t. Recall that the uncertainty in q is at 
worst a factor of 2 while the luminosities differ by a factor of 10 — 
100. The shallow-cusp galaxies thus host an unambiguously larger 
luminosity in their central component. Note that more luminous 
galaxies, Mvt J; —21.5 mag, typically have a shallow cusp, while 
fainter galaxies have steep central cusps. 

3) The 3D half-light radius r^E of the central component of 
shallow-cusp galaxies is typically an order of magnitude larger than 
that of steep-cusp galaxies (Fig|23b. The same is also true for the 
scale radius r_2 (listed in the keys of Figs lll2l"t . which character- 
izes a radius inside which the logarithmic slope of the 3D profile 
falls below —2, and that of the SB profile falls below —1. 

4) The shape parameter n of the inner component of shallow- 
cusp galaxies is always significantly smaller than that of their outer 
component, with the difference reduced for the steep-cusp galaxies 
(Fig|24t. Exceptions, with ricentrai ~ nouter, are NGC4467 and 
NGC4459 whose central region is affected by a huge dust disk, 
and NGC4387 that has ricentrai > nouter, whose light could be 
affected by NGC4406 or, as discussed in section WA] could be a 
3-component system being modelled with 2-components. 

5) All of the nine massive shallow-cusp galaxies 
(AfvT^— 21.71) have n < I for their central component 
(NGC4382 has n=1.55), while for the fourteen steep-cusp galaxies 
usually n > 2.0 and with a larger dispersion; except NGC4478 and 
NGC4434 where n ~ 1 and NGC4551, NGC4458 and VCC1440 
with n ~ 1.5. Lower n implies a larger concentration. Hence, the 
central components of shallow-cusp galaxies appear to be more 
concentrated than those of the steep-cusp galaxies. 

The above observations indicate a new trend with regards to 
the central components of galaxies: Even though steep-cusp galax- 
ies have a higher central density, the central component of the 
shallow-cusp galaxies is far more luminous and massive, spatially 
more extended, and more concentrated than that of the steep-cusp 
galaxies. 

6) The large incidence of n ~ 1 (61 per cent) for the central 
component and especially in all the massive shallow-cusp galaxies 
is interesting. Disks are usually fit with an exponential, or m=l 
Sersic profile. The central regions of these galaxies thus possibly 
support disk-like systems. While this conclusion is very tentative, 
because it is based solely on a SB analysis, in section [84] we de- 
scribe two cases, M87 and NGC4261, whose central regions have 
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Figure 22. Luminosity of components resulting from our multi-component 
DW models as a function of total luminosity of the galaxy (Myy). Galax- 
ies whose component parameters are uncertain (see section l6!2t are marked 
as exceptions in the figure key. The central component of massive (gen- 
erally shallow-cusp) ellipticals are ~ 10^ more luminous than that of the 
less massive (generally steep-cusp) ellipticals. This effect is more apparent 
for the total luminosity in the central-Hintermediate component. Note that 
luminosities of the central and intermediate DW-components are excess lu- 
minosities in those regions with respect to contributions from the outer DW- 
component. A colour version of this figure is available in the online edition. 
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Figure 23. The intrinsic (3D) half-light radii te of all components appear 
to be generally increasing with galaxy luminosity. The r^E of the massive 
shallow cusps are generally a factor of 10 larger than that of the steep-cusp 
galaxies. This trend is stronger for the outer component. Note that, as shown 
in Tabl^ tseo of the outer component is generally slightly larger than r^E 
of the galaxy which is computed using the total hght from all components. 
Also refer to caption of Fig. [22] A colour version of this figure is available 
in the online edition. 
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Figure 24. Einasto index n of components. All massive shallow-cusp galax- 
ies have n<l for their central components, while the steep cusps typically 
have n > 2. The n of the outer component shows a distinct trend increas- 
ing with the total luminosity of the galaxy. Also note the large difference 
in n between the central and outer components of the shallow-cusp galax- 
ies, which diminishes with decreasing luminosity for the steep-cusp family. 
Refer to caption of Fig.|22]for exceptions. A colour version of this figure is 
available in the online edition. 



could recover the total light in the true central component of steep- 
cusp galaxies by fixing, without fitting, the Sersic index of the in- 
ner component to m=\. However, whenever they fit for the Sersic 
index m in the steep-cusp galaxies observed with the HST, they 
typically obtain a wide range of 0.6<m<5.75 for the central com- 
ponent and rarely m~l (refer to the online version of table 3 in 
Hopkins et al. (2009a)). 

This is consistent with the large dispersion in the inner n we 
get for our 2-component DW models in cases where 3-component 
DW models produce significantly better fits. Using three compo- 
nents in these cases, not only is the variance in the inner n reduced, 
we also see a tendency of n~l similar to what Hopkins et.al. find 
for the fits to their simulations with a Sersic profile. 

We caution that only within a limited range of projected 
radii, systems with Einasto profiles of index n~ 1 can be modelled 
with Sersic profiles of index m~l (from equations 3.2 and 3.3 in 
DWIO). Over large radial ranges the Sersic profile fails to model 
a projected Einasto profile. The above comparison should thus be 
taken as qualitative. Nevertheless, this indicates that it may not al- 
ways be meaningful to fit all galaxies with a predecided number 
of components which can further lead to misleading interpretations 
of other structural parameters. It is more meaningful to infer the 
existence of these components either from the residual patterns in 
the data (as in this paper) or motivate their existence through other 
(possibly spectroscopic) observations. 



spectroscopically identified disks, and our 3-component - but not 
2-component - fits show that they have n~l. We reiterate that 
while disks usually have m=l, not all m=l systems should be 
called disks (and vice-versa) without spectroscopic verification. 

To compare with simulations we reviewed Hopkins et al. 
(2008, 2009a,b). While their simulations do not probe the very 
central regions, < 50 pc, they found a similar trend as mentioned 
above; all steep-cusp and shallow-cusp galaxies show a light ex- 
cess in their central regions with respect to the inner extrapolation 
of an outer Sersic profile. They also note that, on an average, they 



8.2 Intermediate component 

Our multi-component models indicate that 14 of the 23 galaxies 
have an intermediate DW-component. These include both steep 
cusp and shallow cusp galaxies. Only in two out of the 14 galax- 
ies (NGC4621 and NGC4434) it identifies features within the outer 
DW-component, but in the rest it is located between the inner and 
outer components, and thus forms a transition region. From Fig|22l 
we also observe that within a factor of two, all galaxies contain 
a similar fraction of the total light in their central-l-intermediate 
DW-components. There appears to be an indication that the mas- 
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sive shallow-cusp galaxies may contain a larger fraction, although 
this is not very clear from our small sample. 

Without additional kinematic or spectroscopic data it is not 
possible to ascertain the physical origin of the intermediate compo- 
nent, but based on the results of existing galaxy formation models, 
we speculate that at least three scenarios are possible: 

We suggest that, the stars making up the intermediate compo- 
n ent - (i) may have formed a s a result of local sta r formation (as 
in lMihos & Hemguistl ( ll994 ). lHopkins et alj j2008l) for the central 
regions), (ii) scattered into this regio n during mergers by a cen- 
tral su p ermassive black hole (SMBH) i Begelman. Blandford. ReesI 



( Il980h . iBovlan-Kolchin. Ma & Ouataerll ( |2004|) ) . or (iii) in the 
case of shallow-cusp galaxies only, stars could have been scattered 
to these radii by the central SMBH, or coalescing binary SMBHs 
whi ch are believed to scour out few 100 pc regions in galaxy cen- 
tres dOualandris & Merrittll2008l) . 

The physical interpretation of the central and outer regions of 
a galaxy depends on how the overall SB profile is modelled. In our 
models for shallow-cusp ('core') galaxies, we observe an excess lu- 
minosity due to the central DW-component with respect to an outer 
DW-component. In section [84l we show that for three 'core' galax- 
ies, the central DW-component does correspond to spectroscopi- 
cally identified real systems. Hence, the central DW-component is 
not necessarily a mere mathematical construct. Alternatively, as has 
been suggested in the literature, the mechanisms mentioned in the 
previous paragaraph are believed to have caused a deficit of mass 
in the central regions, with respect to an inward extrapolation of 
an outer Sersic profile which is fitted to the SB at radii beyond the 
transition, or break radius. In other words, the existence of a shal- 
low cusp may imply - (a) a deficit of mass in the central regions or 
(b) an excess of mass with small n, which is superimposed on to 
the outer DW-component, as in this paper The dynamical interpre- 
tation of the central region will probably be different depending on 
whether (a) or (b) are mostly correct, but either case is consistent 
with SMBH scouring out mass. However, SB analysis alone cannot 
resolve this issue; more data and dynamical modelling are required. 
Further discussion on mass deficits is provided in section [TTI 

8.3 Outer component 

The multi-component DW models reveal a huge dominant outer 
component that usually contains a much larger fraction of the total 
light than any of the other components. Fig. |24] indicates that the 
outer n increases with luminosity. From the SB profiles (Fig. [T]- 
121b we note that for galaxies with outer n>5, the outer component 
makes a significant contribution to the density in the central region. 
This is usually the case for the massive shallow-cusp galaxies, but 
is also seen in smaller, less luminous galaxies, like NGC4473 and 
NGC4434, where the SB appears to indicate a shallow inner slope. 
On the other hand, the outer components of the steep-cusp galaxies 
do not contribute fractionally as much light to the central regions. 

From the above discussion of components in sections 18.118.31 
we observe that all shallow-cusp galaxies have a sharp transition 
to a shallow inner slope. These SB profiles are well modelled by a 
combination of - (i) a central n~ 1 DW-component whose intrinsic 
half-light radius rsE (Table |2ll is greater than about ten times that 
of steep-cusp galaxies, and (ii) a non-negligible contribution to the 
density in the central region from a n>5 outer DW-component. 
The central slope flattening is especially pronounced because an 
Einasto profile with n~l is quite concentrated, compared to that 
of larger n's, and when projected it produces a step-function like 
sharper transition in a log-log plot of density versus radius. 



8.4 Components and systems: Three examples 

The components obtained through fitting ad-hoc functions may not 
correspond to physically distinct kinematic systems or stellar pop- 
ulations, unless the form of our fitting function happens to be cor- 
rect. While a detailed analysis of components and systems for every 
galaxy is beyond the scope of this paper, we present three interest- 
ing connections between the components we deduce from fitting 
the V-band SB profiles , and spec troscopically identified systems. 

Ijaffe et all h993h and iFord et al.l ( 1 19941) reported the earli- 
est detections of nuclear disks around supermassive black holes 
(SMBHs) in the centres of galaxies. We discuss structural simi- 
larities between such spectroscopic detections of systems and the 
components deduced from our multi-component DW models for 
three galaxies in our sample. 

1) Images of M87 (NGC4486) show a prominent nuclear disk. 
The central DW-component of our 3-component model (Fig|2l( 
from the broadband I-band (WFPC F785LP) images, suitably 
scaled to V-magnitudes in KFCB09, has a best fit Einasto shape 
parameter n=1.09 which is consistent with it being a disk. Fur- 
ther, at the scale length of r2 =0.266 arcsec the intensity of the cen- 
tral DW-component drops to 10 per cent; a size c o nsistent with the 
spectroscopic observations o f lHarmsetal1 l ll994h .l Tsvetanov et al.l 
1 1999) study the morphology of the disk in detail with the narrow 
band F658N filter and observe that a significant light excess is de- 
tected inside 0.5 arcsec. We note that at 0.5 arcsec the intensity of 
the central DW-component falls to 1 per cent of maximum. 

2) NGC 4261 also h as large spectroscopically confirmed cen- 
tral systems jjaffe et al.l II993). Ferrarese, Ford & Jaffe (1996)). 
FFJ96 fit a double exponential model with scale lengths of 1.83 
and 8.73 arcsec. Our best fit 3-component model (Fig|6]l show that 
the central and intermediate components have nearly exponential 
profiles. The central component has a shape parameter n=0.9 and 
r_2=2.71 arcsec, while the intermediate component has n= 1.488 
and r_2=9.548 arcsec. 

3) Ano t her w ell studied galaxy is NGC4473 (see 
IPinknev et alj | |2003|) and referenc es therein). SAURON integral- 
field spectroscopic observations temsellem et"ai] \2Q04) shows 
an inner spheroidal system extending through 20 arcsec. Our 
2-component DW model (Fig|14t has a central component with 
n=2.3 and r_2=3.16 arcsec which at 20 arcsec contributes ~ 10 
per cent (2.5 magnitudes fainter) to the total SB at that radius. 

This is a fairly remarkable identification of physically distinct 
systems from a purely surface brightness analysis that had no in- 
formation or priors about the spectroscopic properties of these re- 
gions. While identification of physical systems must not be made 
solely from identification of fit-components through a purely sur- 
face brightness analysis, the above examples indicate that a good fit 
consistent with measurement errors are unlikely to detect spurious 
artificial components. Further, all three galaxies exhibit a shallow 
central cusp. Our models thus allows one to account for structure 
within the shallow-cusp ('core'). This is not possible with either 
the Core-Sersic or Nuker profiles which cannot model deviations 
of the SB profiles from pure power-laws within the central region. 



9 Uniqueness of deprojection 

Given that a multi-component DW-function fits the SB profiles ex- 
tremely well, our next major goal is to explore the conditions under 
which the 3D intrinsic luminosity density profiles can be described 
with a multi-component Einasto model. An attractive feature of 
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such an interpretation is that the Einasto profile can be a likely 
descriptor of both the baryonic, and ACDM N-body dark matter 
haloes, revealing an universality in their functional form. 

Deprojections, however, are generally not unique. In this sec- 
tion we review some of the limitations in obtaining unique depro- 
jections of surface density profiles that have been taken into consid- 
eration while providing the intrinsic 3D luminosity density profiles 
in the next section. 

9.1 Konus and semi-konus densities 



iRvbickil i fT987h showed that based on the angle of inclination with 
the line of sight i of an axisymmetric system, there exists a cone of 
ignorance (0=90° —i), such that a fa mily of densities - called konus 
densities ^Gerhard & Binnevll 19961) — that is non-zero only within 
9, can project to yield zero SB unless the system is seen edge-on 
(i=90°). The range of possible deprojections increases with de- 
creasing i. For triaxial systems this n on-uniqueness increases dra- 
matically jGerhard & Binnevlll996l) . iKochanek & Rvbickil l ll996h 
further extend the range of possible functions through semi-konus 
densities and their linear combinations. 



9.2 Magnitude of effect of semi-konus densities 

Ivan den BoschI l l 19971) (hereafter, vdB97) highlight physical admis- 
sibility conditions that limit the range of possible semi-konus den- 
sities. Through a generalized set of semi-konus densities, he inves- 
tigated how they may (or may not) have an effect on the dynamical 
and photometric properties of galaxies. We summarize here some 
of the key results from his work: 

1) For intrinsic axis ratio, q > 0.8 the effects of semi-konus 
densities are negligible for almost all inclination angles. This is 
similar to saying that spherically symmetric systems have unique 
deprojections. 

2) For oblate spheroids with a constant core, the presence 
of konus densities manifests as detectable wiggles along the mi- 
nor axis, but not along the major axis, an effect seldom seen in 
real galaxies. Although real galaxies seldom have perfect constant 
density cores, this can be an important distinguishing feature for 

shallow-cusp galaxies. 

3) Using a more representative, double power-law dOian et al] 

1 19951) parametrization of the intrinsic 3D density of the central re- 
gions of real galaxies, vdB97 showed that even for an intrinsic axis 
ratio of q = 0.5, the maximum amount of semi-konus density that 
can be added is negligible for inclination angles i > 70°. This is 
seen for a wide range of cusp steepness, —2 ^ a ^ inside the 
core radius, or break radius of the double power-law, say Rt . For 
smaller inclination angles, i < 70°, the maximum semi-konus den- 
sity that can be added increases as the ratio of the scale-length of 
the konus density (say rt) to the core radius of the galaxy decreases 
(refer to his fig. 8). 

A consequence of this effect of decreasing / Rt is that the 
semi-konus density does not add significantly to the mass, or light. 
Furthermore, as r — > the mass of the central supermassive black 
hole will completely overwhelm any contribution of the konus den- 
sity. It should also be noted that the amplitude shown in his fig. 8 
is the ratio of maximum konus-density at r = to the power-law 
galaxy density at r = Rt- Since galaxy densities continue to rise 
for r < Rb, this means the relative strength of konus to galaxy 
density is even smaller at comparable r. 

4) The generalized konus densities are not power-laws. How- 
ever, to characterize their effect in terms of the slopes of the konus 



densities (not to be confused with slope of the galaxy SB profile), 
vdB97 approximates them as power-laws. He showed that if one 
can approximate Pkonus{r) oc r~", then a must be ^ 1. i.e., the 
konus densities themselves can not be too cuspy. 

9.3 The effect of triaxiality 

The discussion in the previous section applies to axisymmetric sys- 
tems. For triaxial systems the range of possible deprojections in- 
creases. However, this increase in non-uniqueness will also depend 
on the degree of triaxiality. 

Triaxiality as well as axisy mmetry can leave an im print on 
the kinematic structure of galaxies. lEmsellem et al J ( [20041 ) provides 
detailed kinematic maps of 4 8 E/SO galaxies from the SAURON 
survey Ide Zeeuw et al.ll2002h ; eight of these galaxies are part of 
our sample. Using an estimate of their specific angular momen- 
tum Xr, within one effective radius ii_E, Jimsellem et al. (20C)3) 
characterized ga l axies as fast (Ah>0.1) or slow (Aii<0.1) rotators. 
ICappellari etal] | |2007|) observe that the fast rotators are typically 
found to be oblate axisymmetric systems and have steep-cusps 
while the slow rotators generally have shallow-cusps and are tri- 
axial. 

Further, lOerhard & BinnevI l ll995h . iMerritt & Fridmatj 



( Il996l) . Iwiuri & MerrittI l ll998l) have shown that triaxial models 
with steep cusps or with a moderately sized supermassive black- 
hole {Mbh ~ 0.005 Algal) are not able to sustain a triaxial shape. 
This indicates that the steep-cusp galaxies which are typically fast 
rotators can be very well approximated as axisymmetric systems. 

For galaxies with shallow cusps which typically rotate 
slowly, theoretically, triaxiality can not be excluded. However, 
ICappellari et al.l ( |2006|) observe that significant triaxiality will lead 
to an increase in the otherwise relatively small scatter in the M/L- 
a relation of the SAURON sample deduced using axisymmetric 
models. It thus appears that the shallow cusps may not be strongly 
triaxial. In our sample the shallow-cusp galaxies also appear nearly 
spherically symmetric in projection and hence we model them as 
oblate axisymmetric systems as well. 



9.4 Non-parametric vs. Parametric deprojection 

In addition to the above iss ues pertaining t o kon us densities, 
iMerritt & Tremblavl i ll 9941) and lGebhardt et alj i ll 9961) highlight the 
importance of non-parametric deprojections. This should be the ap- 
proach of choice to reveal the range of possibilities in 3D since 
small deviations from fits to the SB with ad-hoc fitting functions 
can translate into larger deviations in the 3D distribution. 

However, with very high precision data as used in this paper, 
large deviations, or features in 3D will manifest themselves as de- 
tectable, but possibly small features in the 2D SB profiles. Upon 
projection, smaller features in 3D may remain hidden within mea- 
surement errors. Hence if the goal is to extract smaller, local fea- 
tures in the 3D distribution, then a non-parametric inversion must 
be performed. 

Our goal in this paper is, however, to investigate whether the 
gross properties of the 2D and 3D distribution of light can be de- 
scribed with a multi-component Einasto model. Further, in section 
110. II we show that for two galaxies in our sample that have non- 
parametric deprojections in the literature, our parametric estimates 
are in good agreement with the non-parametric estimates. We hence 
do not take a non-parametric option and note that our inferred in- 
trinsic profile is likely to fall within the confidence intervals of a 
non-parametric deprojection. 
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10 Intrinsic 3D luminosity density profiles 

In this section we discuss the luminosity density profiles of 14 
galaxies in our sample where non-uniqueness of deprojection can 
be reasonably minimized. Based on the discussion in section [93] 
we model all galaxies under the assumption of oblate axisymmetry. 

For oblate axisymmetric systems with intrinsic axis-ratio q 
and with a minor axis inclined at an angle i with the line of sight 
(i = 90° is edge-on), it can be shown that the surface brightness is 
given by: 
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(10.1) 



where, i? is a coordinate along the projected major axis, which for 
oblate axisymmetric systems is the same as the true major axis; 
C, = is the reduced coordinate along the line of sight and r = 

It can be further shown that the observed axis-ratio, g', is given 

by: 



q — \J cos^(i) + sin^(i) 



(10.2) 



When the SB is described in terms of the semi-major axis 
using a spherically symmetric function, generally from equation 
dlO.lt . it can be seen that the inferred intrinsic central density 
po oc Eo(g'/g). An estimate of the light enclosed within an intrin- 
sic radius r, measured along the semi-major axis, when expressed 
in terms of Eo, does not require knowledge of the true axis ratio q 
and inclination angle i, and is given by L{r) oc Eog'. 

For the Einasto profile, equation M.2\ , using ( 15. It and l |5.2t 
this leads to: 
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The p{r) and L{r) profiles in Figures I25l38l have been obtained 
using equation GaUl in Ol and i fTol ). 

Table|2]lists the inclination angles (and the corresponding ref- 
erences) used to obtain the intrinsic luminosity density profiles 
and the cumulative light enclosed as a function of intrinsic 3D 
radius r in Figures 1 251381 For galaxies with observed axis ratio 
> 0.8 for which we could not find inclination angles in the lit- 
erature, we infer their 3D intrinsic luminosity density by assum- 
ing an arbitrary inclination angle i = 85°, to distinguish them 
from galaxies with i = 90° estimated from modelling. Given 
that the residuals of our fits to the 2D surface brightness (median 
rms = 0.032 mag arcsec~'^), are similar to the rms of random 
enws of the data ~ 0.03 mag arcsec^^, we can be fairly con- 
fident that our parametric description of the 3D distribution as a 
multi-component Einasto profile will fall within the wider confi- 
dence limits of non-parametric inversion. This is especially true at 
large R and for galaxies with uniformly low residuals over the en- 
tire dynamic radial range. This is however not true for cases like 
NGC4382 where our 3-component model has systematic residuals 
significantly larger than measurement errors. We also do not pro- 
vide luminosity density profiles for NGC4406 which has a strong 
central dip in its SB profile, as well as galaxies with large inclina- 
tion angles, for example, NGC4261. 



Table 2. Luminosity density sample 



Name 


9^1 IV- 




i 

(deg) 


NGC4472 


0.806 


0.83 


90 


NGC4486 


0.722 


0.96 


90 


NGC4649 


0.828 


0.90 


90 


NGC4365 


0.717 


0.75 


68 


NGC4636 


0.760 




85 


NGC4552 


0.873 


0.96 


90 


NGC4621 


0.742 


0.66 


90 


NGC4478 


0.822 




85 


NGC4434 


0.928 




85 


NGC4473 


0.607 


0.61 


73 


NGC4458 


0.879 


0.88 


90 


NGC4467 


0.813 




85 


VCC1627 


0.928 




85 


VCC1199 


0.869 




85 



Reference 



1 der Marel. Binnev & Davies (1990) 
Cappellari et al. f20071 
Shen & Gebhardt (2010) 
van den Bosch et al. (2008) 

Arbitrary assumed i 
CapDellari et al. (2007) 
Cappellari et al. (2007) 

Arbitrary assumed i 

Arbitrary assumed i 
Cappellari et al. (2007) 
Cappellari et al. (2007) 

Arbitrary assumed i 

Arbitrary assumed i 

Arbitrary assumed i 



NOTES. — Sample of 14 galaxies wiiose intrinsic (3D) luminosity profiles are shown 
in Fias l25l38] Refer to sections |9] and ^| In the columns above: Galaxy type are 
from KFCB09; g^^V" characteristic axis ratio deduced using equation <6.3t : 

q'^^^ and i are the axis-ratio and inchnation angle of the minor-axis with respect to 
the line-of-sight, from the reference listed in the last column. For some galaxies that 
are fairly round but for which we could not find an inclination angle in the literature, 
we have assumed an arbitrary inclination angle of i = 85°. 



10.1 Comparison with non-parametric deprojections 

Some of the galaxies for which we present the intrinsic lumi- 
nosity density, were also deprojected by other authors. The 2- 
component model for NGC4649, Fig. |27] can be compared to the 
non-para meteric inversion of the sa me data set (from KFCB09) pre- 
sented in IShen & GebhardJ feOlOl) (SGIO). Note that in this pa- 
per we use values for observed axis ratio g and inclination angle i 
quoted by SGIO. The overall profiles are consistent once a correc- 
tion is made for the distance to the galaxy; SGIO adopt 15.7 Mpc, 
while we use 17.3 Mpc. The central luminosity density is steeper in 
SGIO. However, one must note that the central-most points in the 
surface brightness have larger errors, arising primarily from psf- 
deconvolution. 

Also as mentioned in section (6] including a third component 
improves the fit to the SB in the central regions (Fig.O, making it 
consistent with the ~ 0.1 mag arcsec~^ errors in the central-most 
data points. We do not use the 3-component model because the F- 
test rejects it at 1.2a. One can thus use the 3-component model 
purely as a fitting function (without making strong inferences on 
the resulting components) and in doing so we find very good agree- 
ment with the non-parametric deprojection in SGIO. This indicates 
that a subdued third component could well exist, as shown in Fig.[3] 

A similar comparison of our 3-component parametric model 
(Fig.[2l> with a non-paramet ric deprojection of NGC4486 (M87) in 
iGebhardt & Thomas! l [2009l) also shows that they are consistent. 



10.2 The form of the 3D intrinsic structure 

Most of the galaxies for which we present luminosity density pro- 
files are fairly round: ten are E0-E2, and four are E3/E4. Twelve 
galaxies have large, ~ 90° inclination angles, and two are at i 
~ 70°. None have any wiggles or dips in the central regions of 
their SB profiles. These considerations ensure that the contribu- 
tion of semi-konus densities to the density profile is small, or non- 
existent. Further, our multi-component models have residuals con- 
sistent with measurement errors over a large dynamic radial range 
that allows for a parametric deprojection. 
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Figure 25. Multi-component Einasto models of the intrinsic (3D) volume luminosity density profiles for 14 galaxies in our sample (see Table|2]and 
section [Tol assuming oblate-axisymmetry. Only the statistically significant best-fitting models (Table[T) are shown. The component profiles have been 
computed using equation 110.31 in \Y.2\ . Colours and line types are as in Fig[T] The figure keys list the Einasto shape parameter n and the intrinsic (3D) 
effective or half-light radius of the components (in kpc), estimated from the best-fitting values of r_2, using equation lS.SI The total half-light radii in 
kpc, intrinsic ir-jE) and projected (R2E), are also shown separately in the figure panel. Also listed are the observed axis ratio q', the inchnation angle 
i of the minor axis to the line-of-sigh t wit h suffixes labelling the re ferences (see Table[2) - DW (this p aper), vM (van der Mai'el et al. 1990), SN stands 
for SAURON jCappellari et al.ll2007h . SO fehen & GebhardllEoiol) and vB jvan den Bosch et alj2008l) . For some galaxies that are fairly round and for 
which we could not find an inclination angle in the literature, we have assumed an arbitrary i = 85°. These cases are labelled as iarb- Generally we 
use q' from the same reference that contains the i listed in Table|2] However, if <7^^y ~ q' of the reference, we use q^j^. The intrinsic axis ratio qint 
is computed using equation 1 10.21 The horizontal axis showing the intrinsic (3D) radius (in kpc) is up to ~ 1.5 X the projected radius of available data. 
Above NGC4472 (left) and NGC4486 (right). (Colour versions of these figures are available in the online edition.) 



We conclude that for the 14 galaxies presented in this section, the 3D intrinsic density profile can be described with a multi-component 
Einasto model. These galaxies span a wide range of luminosities —24 < Mvt < —15 and belong to both the steep-cusp and 
shallow-cusp families. It is therefore likely that the intrinsic 3D baryonic density of other ellipticals can also be described with a 
multi-component Einasto model. 



10.2.1 Considerations for detailed modelling 

Here we note a few important factors that may affect our interpretation of the intrinsic luminosity profiles. 

i) We have assumed a constant ellipticity while galaxies seldom have constant ellipticity. A varying ellipticity can be incorporated in 
our multi-component models for a more accurate deprojection. 

ii) We show luminosity density profiles of 14 galaxies which we could justify as oblate axisymmetric systems. However, some galaxies 
are prolate while some are triaxial and some have a combination of axisymmetric and triaxial regions. For such galaxies, the 3D 
intrinsic density should not be inferred from a SB analysis alone. One may assume a multi-component Einasto profile for the intrinsic 
density and use additional kinematic information to include triaxiality in the kinematics to constrain the Einasto profile parameters 
either in 3D or in 2D through the DW-function. 

iii) The luminosity density profiles are subject to uncertainties in estimating the absolute magnitude (which are subject to uncertainties 
in distance measurements) and Galactic extinction. These uncertainties however do not affect the overall shape of the density profile. 
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Figure 26. Cumulative light enclosed within an intrinsic (3D) radius r{kpc) for the 14 galaxies for which we infer the intrinsic luminosity density (see Table|2] 
and section [Tot . The component contributions have been estimated using IIP. 4) . Profiles are shown only for the statistically significant best-fitting models 
(Table[T). The observed axis ratio q' and the total half-light radii, intrinsic r^E ™d projected R2E, sre also shown separately within the figure panels. Colours 
and line types are as in Fig[T]and details on figure keys and labels are as in Fig |25l Above NGC4472 (left) and NGC4486 (right). (Colour versions of these 
figures are available in the onhne edition.) 
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Figure 27. Intrinsic (3D) luminosity density for NGC4649 (left) and NGC4365 (right). Refer to caption of Fig[25]for details. 
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Figure 29. Intrinsic (3D) luminosity density for NGC4636 (left) and NGC4552 (right). The vertical marker on the profile for NGC4552 shows the inner radius 
beyond which the luininosity density profile shown, can be trusted. This is because the best fitting 3-component DW model for this galaxy (Fig llO) is not a 
good representati on of the surface brig h tness profile within 0.0 45 arcsec where the light is affected by the variable UV flare activity inteipreted to arise from 
a low-level AON iRenzini et aP il995r) . ICapDeflari et all h999h ). Also refer to caption of Fig|25]for details. 
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Figure 30. Cumulative volume luminosity profile for NGC4636(left) and NGC4552(right). The vertical marker on the profile for NGC4552 shows the in- 
ner radius beyond which the profile shown, can be trusted. This is because the best fitting 3-component DW model for this galaxy (Fig llOt is not a good 
representation of the surface brightness profile within 0.045 arcsec. Refer to caption of Fig [26l and Fig|29]for details. 




Figure 31. Infi-insic (3D) luminosity density for NGC4621 (left) and NGC4478 (right). Refer to caption of Fig|25]for details. 
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Figure 32. Cumulative volume luminosity profile for NGC4621 (left) and NGC4478 (right). Refer to caption of Fig[26]for details. 
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Figure 33. Intrinsic (3D) luminosity density for NGC4434 (left) and NGC4473 (right). Refer to caption of Fig|25]for details. 
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Figure 34. Cumulative volume luminosity profile for NGC4434 (left) and NGC4473 (right). Refer to caption of Fig[26]for details. 
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Figure 35. Intrinsic (3D) luminosity density for NGC4458 (left) and NGC4467 (right). Refer to caption of Fig|25]for details. 
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Figure 36. Cumulative volume luminosity profile for NGC4458 (left) and NGC4467 (right). Refer to caption of Fig[26]for details. 
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Figure 37. Intrinsic (3D) luminosity density for VCC1627 (left) and VCCl 199 (right). Refer to caption of Fig|25]for details. 



34 Barun Kumar Dhar and Liliya L.R. Williams 




11 Mass deficit in massive ellipticals 

IPeeblej [WT^ predicted that a compact massive object, of mass 
M, can alter tiie density structure around its radius of influence, 

ri„f=GM/al (11.1) 

such that for r < ri„f a stellar density cusp (in 3D) cx r~^/* 
forms with a velocity dispersion a{r) oc r~^^^ ; and a core, of con- 
stant density and velocity dispersion, forms within a core-radius 
rc »rinf . Here (j« is the veloci ty dispersion at r»rinf. 

iBegelman. Blandford. Reei ( fUsQ) discuss the formation, 
evolution, eventual coalesence and possible recoil, of a binary 
SMBH, as a result of galaxy mergers. Since then, N-body simu- 
lations have shown that the evolution of the binary leads to a co- 
evolution of the density profile of the galaxy, in a region around the 
radius of influence, due to ejection of an amount of mass propor- 
tional to the mass of the binary. 

Our work, describing the structure of ellipticals as a superpo- 
sition of DW-profiles, implies the presence of excess light in the 
centre of all ellipticals with respect to an inner extrapolation of the 
outer DW-components. We emphasize that our models describe the 
current structure, post any modification by the binary, and can be 
consistent with mass ejection by SMBHs. 

In this section we discuss a related concept of mass deficit in 
massive ellipticals, which has its genesis in - 

i) the identification of massive ellipticals exhibiting a shal- 
lower surface brightness (SB) profile in their central regions, in- 
wards of a projected break radius, with respect to the trend of SB 
outside the break radius - i.e., a light deficit. This is in contrast to 
less massive ellipticals that seem to contain steeper light profiles 
inside the break radius - i.e., a light excess; 

ii) the assumption that evolution of the binary SMBH is the 
sole factor responsible for the observed break radius; and 



iii) the assumption that in the absence of the binary SMBH, 
the trend (functional form) of SB outside the break radius would 
have continued its form all the way to the centre of the galaxy. 

In section 17 1.3 1 we show that estimates of 'observed' deficit 
in real galaxies, made with respect to an inner extrapolation of the 
SB profile outside the break radius, not only have large variation 
with respect to predictions from N-body simulations but also vary 
widely between different researchers for the same galaxy. We also 
compare signatures of mass-ejection by the binary SMBH with the 
observed profiles of real massive ellipticals and find that current 
predictions from N-body simulations are not able to account for 
some of the largest 'cores' in Virgo galaxies. Prior to such discus- 
sion, in section fTl.il we briefly describe the various phases of evo- 
lution of the binary SMBH and in section [Tl.2l we review the pre- 
dictions from simulations about the amount of mass ejected, spa- 
tial extent and resulting slope of density profiles, in the context of 
'dry' dissipationless (gas-free) mergers, believed to be the forma- 
tion mechanism for the massive 'core' ellipticals. The discussion is 
thus directed towards systems with the most massive > 10* Mq 
SMBHs, as is the case for the 'core' ellipticals in Virgo. 



11.1 Phases of Evolution 

The significant phases in the evolution of the binary SMBH are: 

Phase 1: Subsequent to a galaxy-galaxy merger, a SMBH bi- 
nary is said to form when the SMBH separation reduces to a dis- 
tance Tinf - the radius of influence of the more massive SMBH. 

Phase 2: Their separation shrinks as the binary loses energy 
due to dynamical friction and by scattering of stars through the 
gravitational sling-shot mechanism. A 'hard-binary' is eventually 
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defined to form when the separation reduces to a size, 



^hard 



GfJ. ^ q 
Aal 4(l + q)2 



(11.2) 



where, jjh =M\M2/M\2 is the reduced mass of the binary with 
q =M2 /Ml ^ 1 and Mi2=Mi + M2 is the combined total mass of 
the binary, ahard is thus a fraction of ri„/ (equation (TTTTJ) 

At this stage, due to depletion of stars, the decrease in 
binary separation can stall unless additional mechanisms can 
continue remov ing energy f r om th e binary to drive it to - 
wards Phase 3 . iBercziketalJ j2006l). iHoffman & Loebl j2007h . 
iBerentzen et al.l ( l2009h . iKhan. Just & Merritt I j201lh show mech- 
anisms through which s talling may easily be avoided. Also, see 
iDotti. Sesana & Decarlil | |201 li) for a review and references therein. 

Phase 3: The binary then continues to harden until its semi- 
major axis reduces to a size where energy loss due to emission 
of gravitational waves begins to dominate energy losses due to scat- 
tering of stars, driving the binary towards coalesence. 

Phase 4: During the last few stable orbits before coalesence, 
the centre of mass recoils with a recoil or kick velocity that sen- 
sitively depends on the binary mass-ratio and spins, relative spin 
allignement and orientation of spin-axis with the angular momen- 
tum vector. If the kick velocity Vkick > Vesc, the escape velocity, 
the coalesced SMBH can be completely removed from the galaxy. 
For lower kick velocities, the SMBH performs damped oscillatory 
motion in the central regions and eventually settles down to a Brow- 
nian motion about the centre. 



11.2 Results of Simulations 

11.2.1 Mergers in hierarchical models of galaxy formation 

Volonteri. Haardt & Madaul l l2003h (V-l-03a) and 

Volonteri. Madau & Hardt] l l2003t) (V-l-03b) have shown that, 
in galaxy mergers: 

i) SMBHs of mass ratio q—l occur at high redshift 2^20, 
while for 0<2<14, q<0.2 with typical values of g~0. 1. 

ii) Assuming that the progenitors start with p oc r^^ pro- 
files and mergers (accompanied by coalesence of the binary) not 
only erase the cusps completely to form a constant density core 
but also preserve such cores, they estimate that at the end of a suc- 
cession of mergers the resulting mass deficit, Mdef, is given by 
Md<./=8.2±3.8) Mbh (equation 15 of V-l-03b). Note that most 
'core' galaxies have a shallow cusp rather than a constant density 
core. The above estimate should therefore be considered a limiting 
value. At the same time, mass deficit due to the scouring effect of 
SMBHs in Phase 4 is not accounted for in this study. 

iii) Even when the entire cusp is erased and a constant 
density core forms, the resulting core size is only ~60 pc for a 
1O"M0 halo (fig.l and 2 of V-l-03b) ~ a mass resembling the 
most massive 'core' ellipticals in Virgo. For example, assuming 
Mhaio/ Lv=5Q the halo mass of M87, the giant 'core' galaxy 
in Virgo with Lv-t = 1-53x IO^M^q, is ~7.65x IO^^Mq. The 
break radius for this galaxy is however ~600 pc>>60 pc core 
formed by the merging binaries; a point we will return to in the 
next section. V-l-03b also find that such core sizes are typically 
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11.2.2 N-body estimates of Mde/ before coalesence 

Using iDehned l ll993h models - a power-law of the form p(r) oc 
r ' at small r - to describe the intial density profiles, iMerritJ 



( I2OO6I) conducted 39 N-body simulations to study the evolution 
upto the hard-binary stage (Phase 2), for the case of 'dry' merg- 
ers in non-rotating spherically symmetric systems, including nine 
remergers. He showed that, after each merger: 

i) Mdef~0Jq'^'^Mi2- For initial density profiles as Dehnen 
7=0.5 model and for typical values of g=(0. 1,0.25) (V-l-03a), 
Mde/=(0.33,0.46)A/i2 (table 1, Merritt06). 

ii) there is a lowering of the density profile in a region about 
the size of the radius of influence. This can be seen by comparing 
r/^ from his table 1 with the radius at which reduction in density 
occurs after {he first merger in his figs. 6a,6b and 6c. 

[Merritt. Mikkola & Szell ( 2007), MMS07, study the evolution 
in Phase 3. For intial Dehnen density profiles with 7=0.5, and 
for g=0.1, and Mi2 = 10*Mq, they obtain Mde/=0.9Mi2. When 
added to the mass deficit upto Phase 2, this yields Mde/~1.2Mi2. 
Further, as can be seen from their figure 21, the net deficit at the 
end of Phase 3 is also within a region ~ - the region of influence 
defined as the radius containing a stellar mass equal to twice the 
total black hole mass. Note that rh is typically much less than the 
observed break radius ('core') of massive ellipticals. 



11.2.3 N-body estimates of Mdef after coalesence 

Assuming that the binary coalesces and recoils with a kick 
velocity (Ph a se 4) iBovlan-Kolchin. Ma & Ouataerj ( |2004|) and 
iMerritt et alj j2004h investigate the resulting effects on the den- 
sity structure. They find that for Vfcicfc>(0.25-035)Vesc, cores of 
size Vh can form as a result of the kicks; with larger deficits for 

iGualandris & MerrittI ( l2008h (GM08) explore this phase in de- 
tail using a Core-Sersic (CS) profile (in 2D) and its deprojection 
jTerzic & Grahamll2005l) for the density profiles. They found that: 

i) For 0.3 Vesc ^ Vkick < Vesc, the coalesced SMBH performs 
damped oscillations about the galactic centre and loses energy dur- 
ing each passage through the pre-existing 'core'; the later assumed 
to be the break radius of the CS profile. 

ii) the flattest profiles with CS 7<0.05, leading to the largest 
deficit, occurs for Vkick ^ 0.8 Vesc (their fig. 12). It is important to 
note that such flat profiles are produced given an initial profile with 
a fairly low 7=0.55. i.e., it is not clear whether such flat profiles 
can be produced for larger initial 7, even with near Vesc kicks. 

iii) the resulting Mdef/Mi2 ~ 5.08 {Vkick /Vescy-'^^ , in 
Phase 4, and depends weakly on Mi2/Mgai', where Mgai is the 
mass of the galaxy. Consequently for Vkick = (0.1,0.4,0.8)Veac, 
Mdef = (O.I,I.O,3.44)Mi2. Further, Vkick /Vesc decreases with 
galaxy luminosity. Hence the lower values of Mdef are more likely 
for the massive shallow cusp ('core') galaxies; for example, GM08 
show that Vkick=0-^Vesc=550km s^^ is a likely kick-velocity for 
the most massive core galaxy in Virgo, NGC4472 (M49). 

Self-consistently combining A'Idef from all phases, for merg- 
ers of ~10*A/q SMBHs with typical g=0.1 and Vk^ck/Vesc=OA, 
the total cumulative A-Idef is I.IMbh per merger. We shall show 
in the next section that estimates of Mdef of 10,20 and sometimes 
AQMbh, from high-resolution SB profiles, cannnot be explained 
by the above results of N-body simulations. 

iv) and, the radius of the pre-existing 'core' (here, break ra- 
dius) expands by 15 per cent for Vkick=Q-'2Vesc to about 70 per 
cent for Vkick=Q-9Vesc and is lesser for larger Afs_f/ (refer to their 
table 1 and 3). Note that in the simulations upto Phase 3, the ra- 
dial extent of the core (~rh) formed by the SMBHs is less than the 
break radius of massive 'core' ellipticals. Since the break radius 
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has been assumed to correspond to the core formed by the binary 
through Phase 3, the fractional expansion in size of the break radius 
is more relevant than the absolute eventual size of the break radius. 

11.3 Estimates of observed mass deficit 

In this section we investigate whether the signatures of mass deficit 
in N-body simulations are in congruence with the estimates of mass 
deficit from SB profiles of galaxies and highlight limitations of ex- 
isting models of galaxy structure - a single sersic, a core-sersic and 
a Nuker profile - in estimating such mass deficits as well as in con- 
firming predictions from simulations. 

In real galaxies, estimates of the amount of mass deficit and 
the extent of break radius vary widely for the same gal axy. This 
can be seen by comparing the mass deficits computed in lGrahamI 
( |2004|) (their table 1), with data from Rest et al.(2001). 

For instance, for NGC4168, Graham04 estimated a projected 
break radius, Rb, of 0.72 arcsec (108 pc) and a Mdef=l-2MBH 
using a CS profile, while estimates with a Nuker profile yield 
7?6=2.02 arcsec (3 03 pc) and a Mdef=23 . 5MBH ', for the later, 
the prescription in iMilosavlievic & Merrit assuming that 

the unscoured intrinsic (3D) profile within the break radius was 
p{r)(xr^^, has been used. However, IMilosavlievic et ahl j2002l) 
(MMRvB02) estimate an i?6=1250 pc (their table 1) using the 
same data of Rest et al., and obtain Mdef='^^MBH ■ Not only do 
these estimates differ by a factor of 40 for the same galaxy, even 
the most extreme estimates of Mdej from N-body simulations are 
not able to account for deficits as large as (20-50)Mfl£r . 

For NGC2986, MMRvB02 obtain a lesser estimate of 
Mdef=5.24A4BH from a larger break radius of 400 pc 
compared to the (Nuker,CS) profile estimates (Graham04) of 
Mdef =(26.7,7.02)A/sH over a smaller Rt =(174,94) pc. 

Assuming initial profiles, MMRvB02 show that on av- 
erage i\/de/ = 10 Mbh for the case of ellipticals. This estimate is 
similar to the predictions of V-l-03a (~ 8 Mbh) and the average 
value of 10 Mbh in KFCB09 for Virgo ellipticals . However they 
do not agree with the (2.4±1.8)MBjj estimates of iFerrarese et al.l 
( I2OO6I) (F-l-06) for the Virgo ellipticals, the (2 1± 1.1 )MB_tf es- 
timate of Graham04, the 2 Mbh estimate of iHoffman & Loebl 
( I2OO7I) (HL07), including Virgo elliptic als, and the (2.29 ±0.67) 
Mbh and (1.24±0.3) Mbh estimates of lHvde et all j2008l) . 

For the giant 'core' elliptical in Virgo, M87(NGC4486), 
HL07 estimate 2.5 Mgj/, MMRvB02 estimate 8.7 Mbh and 
iKormendv & Bended | |2009|) (KB09) estimate 14.1 Mbh, all 
scaled to Mbh =3.5 x 10^ Mq. Despite such large deficits, the SB 
profile of M87 shows a distinct rising trend; apart from the fact that 
it has an unsually large break radius. NGC4649 also has a fairly 
large break radius. While KB09 estimate 5Mbh, HL07 estimate 
only 1 Mbh, both for Mbh =2x10^ Mq. Another Virgo galaxy 
with an apparen t 'core' is NGC4382. K B09 estimate a deficit of 13 
Mbh, however I Gultekin et alj j201 ih find that this galaxy has an 
unusually small black hole mass, consistent with no black hole, yet 
obtain a mass deficit of 45.6 Mbh- 

Estimates of Mdef /Mbh do depend on a number of factors 
like uncertainties in estimates of Mbh, the assumed form of M/L 
etc. However, the large differences are generally due to discrep- 
ancy in ~ estimation of Rt; estimation of the inner power-law in- 
dex; and the huge uncertainty in our assumptions about what might 
have been the shape of the density profile prior to the action of the 
SMBHs. For instance, KFCB09 defines the central region of light 
deficit or 'extra-light' based on the region over which a Sersic fits 
their estimate of the non-central region. F06 on the other hand fits 



a Core-Sersic model which defines the central region as the region 
inside Rt of the core-sersic model. Consequently the two meth- 
ods yield widely varying estimates of what the 'central region' of 
deficit is. This leads to a large discrepancy in the estimated Mdef- 
Als o see fig.2 of MMRvB0 2 

iHopkins & Hemguistl | |201(]|) (HH10) suggest a new non- 
parametric method to obtain mass deficit and find that Mdef varies 
with radius such that at ~ 100 pc, Mdef /Mbh ' assymptotes to 
a maximimi of 0.5-2' and for the largest galaxies at core-radius of 
~ kpc, Mdef~(2-4)MBH- While, this is similar to the estimates 
of F-l-06 and Graham04, it is at odds with the KFCB09 estimates 
where the authors reasoned that Mdef ~ 10 Mbh accounts for the 
net mass deficit through Phase 4 (results from MMS07 and GM08) 
and in agreement with V-l-03a and MMRvB02. 

As discussed in section fTl.2.31 Mdef=(lO-20)MBH are hard 
to explain using the results of N-body simulations. However, even 
if such were to be true, the resulting region over which KFCB09 
and MMRvB02 estimates such deficit is much larger than the re- 
gion of influence of the SMBH over which simulations predict mass 
deficits; a point also noted in MMRvB02. This is also true for esti- 
mates of Rt using a Core-Sersic or Nuker model; especially for the 
largest 'cores' as in NGC4486 and NGC4649. The shallow cusps 
('core') are thus not entirely due to mass ejection by the SMBHs. 

11.4 Deficit or excess? 

Disagreements in the literature are limited not only to the amount, 
but also the sign of the deficit. In MMRvB02, NGC4478 is shown 
to have a mass deficit of 15.85 Mbh- However, its SB profile 
shows an apparent 'cusp' (see Fig[T5] of this paper), due to which 
F+06, Cote et al.(2007) (C+07), KFCB09 infer this galaxy to have 
a light (and mass) excess. Similarly in MMRvB02, NGC4473 has 
the highest Mdef=22-'i9MBH amongst the Virgo ellipticals. F-l-06 
and C-F07 also conclude that this galaxy has a mass deficit. But 
KFCB09, argue that this galaxy has 'extra-light' and not a deficit. 
Further, in this paper and in Hopkins et al. (2009 a,b), the authors 
argue that all galaxies can be modelled with extra light above an 
inner extrapolation of the outer components. 

Thus, whether the SB profile of a galaxy 'exhibits' an excess 
or deficit also depends on ones methodology. 

11.5 Discussion 

Given the large disagreements in estimates of Mdef, the region af- 
fected by the binary SMBHs; and which profile best describes the 
light of ellipticals; it can not be said that robust estimates of mass 
deficits, in real galaxies, have been made. Consequently, validat- 
ing predictions of N-body simulations and especially the stage(s) 
of SMBH evolution responsible for the observed deficit turn out to 
be largely uncertain. 

Further, extrapolating the profile from the break radius, Rt, 
presents a number of conceptual difficulties: 

i) When the SB profiles of the massive 'core' galaxies are 
modelled with the Core-Sersic (F-l-06) or single Sersic (KFCB09) 
profiles, the non-central regions, beyond Rt, generally have Ser- 
sic index m > 5; sometimes as large as 9 and 11. Larger the m, 
steeper the density as i? — and consequently larger the extrap- 
olated density. Since estimates of Mdef are correlated with an es- 
timate of Rt and in most cases Rt is much larger than the region 
of influence of the SMBHs (example: NGC4486, NGC4649), ex- 
trapolating profiles with large Sersic indices will invariably imply 
larger, but incorrect, mass deficits (also see HHIO). 
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ii) The SMBHs are believed to affect only the central regions, 
around the radius of influence, and not the global profile of the 
galaxy. While it is believed that the mass deficit is due to ero- 
sion of the steep cusps, it remains to be explained what forms the 
large m~9 Sersic profiles, in the outer regions of only the 'core' 
galaxies. If SMBHs have not influenced these regions, some other 
processes have shaped the formation and evolution of the 'core' 
galaxies. Their shallow cusps could well be influenced by these 
processes as well and not entirelly due to the core-scouring effect. 

Infact, F-l-06, KFCB09 and iLauer et alj ( l2007h observe that 
shallow- and steep- cusp galaxies have different characteristics, not 
just profile shapes, and are likely to have had different formation 
pathways. Their central regions could well have fundamentally dif- 
ferent initial shapes as well, than being a mere extrapolation of 
the profile outside _R(,. Further, baryonic effects are usually more 
dominant (over dark matter) in the central regions of all galaxies - 
core and cuspy alike - leading to relatively more concentrated cen- 
tral, than outer, regions (or components). This is precisely what our 
modelling indicates (section [STI and Fig|24b. If this is true, it will 
not be meaningful to extrapolate the outer component inwards to 
estimate the 'mass deficit' and neither would comparing possibly 
fundamentally different profiles of 'cuspy' and 'core' galaxies, to 
get non-parametric estimates of mass deficits, as in HHIO. 

iii) If the massive ellipticals have formed through multiple 
mergers, it is likely that their outer profile has also been built up as a 
cumulative effect. Consequently their Einasto or Sersic index ought 
to have evolved with their merger histories. Recall that the large in- 
dices, for the massive galaxies, are due to the gradual fall-off of 
their extended stellar light well outside the central regions. We are 
inclined to believe that this index was different in the past and has 
evolved over the merger history of these galaxies. The large spatial 
extent of these galaxies could also be indicative of a larger number 
of mergers, compared to the much smaller steep cusp galaxies; a 
point noted by KFCB09 as well. If this is true, and mergers alter 
the shape (Sersic index) of the outer profile, then extrapolating the 
current observed outer Sersic profile (index) to estimate the initial 
density profile will give us an incorrect estimate of the mass deficit. 

In this paper, we have shown that all galaxies, shallow and 
steep cusps, can be modelled with very high precision over a large 
dynamic range, as a superposition of DW-components and that they 
all have a 'light excess' in their central regions with respect to an in- 
ner extrapolation of their outer components. As shown in sectior[7] 
not only do the DW-models fit the SB profiles better than the Core- 
Sersic and Nuker profiles, at least for three galaxies the central 
component correlates with a real physical system and is not a mere 
mathematical construction of the fitting process. 

Mass ejection and consequently some deficit due to core- 
scouring binary SMBHs could well have occurred and these are 
likely to have shaped the central components of our DW models 
as well. For instance, in this work, we find that the central compo- 
nent of all 'core' galaxies have an Einasto index 7i~I. It will be 
instructive to see if such components can be robustly isolated in 
SPH+N-body simulations and how well can they be fit with n~ 1 
Einasto profiles. Similarly, as noted in section[8j2] it is possible that 
the intermediate DW-components within the central kiloparsec or 
so, could have formed as a result of evolution of the binary SMBHs. 

Recent advances in N-body simulations with large particle 
numbers is a promising development that can shed much light on 
the dynamical role played by the SMBHs. Unlike simulations, how- 
ever, observations of real galaxies do not have the advantage of 
knowing what the unsecured profile was. Hence, estimates of such 



effects from observations should be made with caution and with 
respect to a robust model for galaxy structure. 

From the above discussion, we conclude that mass ejection 
due to SMBHs is unlikely to be the sole cause of the shallow-cusps 
in massive ellipticals. The role of processes shaping the global 
structure of galaxies should also be accounted for. 



12 The outer n of ellipticals 

In this section, we present two empirical speculations about the 
structure and formation of elliptical galaxies. 

(1) Our multi-component models reveal an interesting prop- 
erty of the most luminous, ^^IO^^L^q galaxies in our sample. The 
Einasto shape parameter n of their outer component is very similar 
to the n of pure dark matter haloes; 5 < n < 8, as shown in Fig.|39l 
We remind the reader that because dark matter haloes are well fit 
with Einasto profiles (N-l-04, M-l-06), a direct comparison with the 
n of our galaxies is possible0 

The plus symbols in Fig.|39|represent 26 galaxy-size dark mat- 
ter haloes compiled from the high resolution AC DM N-body simu - 
lations of loie mand, Moore & Stadel (2004), N avarro et al.l(l2004l) . 
Merritt et aTj l l2006h. iPrada et al.l ( |2006|) . IStadel et all ( l2009h and 
Navarro et all l l2010h . We assume that M200 of these haloes is 
80 per cent of Mgau the combined mass of the dark matter and 



baryons. To translate the stellar mass of our Virgo galaxies to Mgai 
we have assumed (i) a stellar mass-light ratio M*/Lv't=5 for all 
galaxies, and (ii) the stellar mass comprises 20 per cent of Mgai- 
Observe that the masses of the galaxies hosted by these dark matter 
haloes, ^IO'^^-IO^^Mq, are comparable to those of the low and 
high luminosity ellipticals in our sample, — 24< MvT< — i5. 

If we interpret the dark matter simulation results to mean 
that Einasto shape parameter 5 < n < 8 characterizes collision- 
less relaxation, then we can tentatively conclude that the distri- 
bution of stars in the outer regions of massive ellipticals, which 
have Uouter^TiN-body, was also shaped by collisionless processes. 
This conclusion is broadly consistent with the prevailing notion that 

massive ell ipticals form through dry, di s sipatio nless mergers. 

(2) iFerreras. Saha & Williain3 ( |2005[) and iFerreras et all 

l l2007h have shown using a combination of strong lensing and stel- 
lar population synthesis models that the outer regions of 
massive galaxie s are more dark matter dominated than those of 
smaller galaxies. ICappellari et al. I( l2006l) arrive at a similar conclu- 
sion that fast-rotating galaxies, which are generally low-luminosity 
with steep cusps, have relatively lower dark matter content than 
the slow-rotating galaxies which ar e usually massive with shal- 
low cusps. Also, lAuger et"al] l l2010h observe that the mean dark 
matter fraction within -Rb/2 increases with galaxy size and mass. 
These observations are consistent with Planetary Nebula observa- 
tions (D ouglas et al. (2007), Napohtano et al. (2007,2009,2011), 
lTortoraetal.l ( l2009t) ) that the outer regions of small and intermedi- 
ate mass ellipticals have varying degree of low dark matter content, 
while dark matter is quite dominant in the outer regions of massive 
ellipticals. These observations indicate that all massive ellipticals 
are dark matter dominated. 



t Papers describing fits to the haloes in N-body simulations use the recip- 
rocal of our shape parameter, a=l/n, and typical a values hover around 
0.17 which corresponds to n=5.88. 
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Figure 39. Einasto index n of ACDM N-body haloes and that of the outer 
component of the Virgo galaxies in this paper. The outer component of 
the most luminous (massive) galaxies have an n similar to that of the N- 
body haloes. This is not true for most of the lower luminosity ellipticals. 
The solid box at n=5.53 is NGC4434 and the open circle at n=9.56 is 
NGC4621 whose outer-n may be uncertain (section [6^ . The other solid 
box at n=4.98 is VCC1440, a low luminosity galaxy bordering the dwaif- 
elliptical population. Not included are NGC4406 and NGC4382, whose 
outer components are severely perturbed. Refer to text ( section [T2t for dis- 
cussion. A colour version of this figure is available in the online edition. 

Since all the massive ellipticals in our sample have outer 5 < 
n < 8, this could indicate that such galaxies are likely to be dark 
matter dominated. There might even be a relation between the outer 
n of the galaxy's SB profile and its dark matter content. 



13 Summary and Discussion 

Critical to our understanding of galaxies, their formation and evo- 
lution, is our ability to accurately quantify the galaxies' structural 
properties and their variation with mass. The Sersic profile is the 
single most commonly used function for describing galaxy light 
distribution. The works of Caon, Capaccioli, Einasto, Ferrarese, 
Graham, Kormendy, Lauer and their collaborators have revolution- 
ized our understanding of galaxy structure. These authors have 
shown us that a single Sersic profile does not fit the surface bright- 
ness distribution of ellipticals consistent with measurement errors 
over a radial range larger than 2-3 decades. To extend the radial 
range of the fit, especially in to the central regions, the core-Sersic 
(6 parameters) and the Nuker (5 parameters for the central 10-20 
arcsec) profiles have been introduced, with a further addition of a 
King model (3 more parameters) for the nuclear region. 

But even with these flexible and multi-parameter models, the 
fit residuals often exceed the measurement errors in some radial 
ranges, in spite of the overall rms of residuals being low due to the 
large regions over which the Sersic profile fits well. Double Sersic 
models (Hopkins et al. 2009 a,b) provide an improvement, but still 
the residuals remain larger than measurement errors. 

In addition to fit residuals exceeding measurement errors, the 
Sersic profile has another drawback. The galaxy structure, dynam- 
ics and evolution exist in 3D, so it is more meaningful to describe 
galaxies using 3D functions. The Sersic profile is an intrinsically 
2D distribution whose deprojections preassumes an infinite 3D ex- 
tent. If, for example, a galaxy has a truncation radius in 3D, a Sersic 
profile will not be able to model that. The parametric forms for de- 



projecting a Sersic profile (and power-laws, for that matter) are not 
well behaved near the centre. Further, these analytic functions are 
either reasonable but not extremely accurate (PS97, LGM99), or 
are accurate but extremely complicated (BGl 1). 

To overcome the above limitations, we propose that the 3D lu- 
minosity density may be described with a multi-component Einasto 
profile whose parameters can be estimated by modelling the surface 
brightness using a multi-component DW- function. This is similar to 
the observations of Einasto and collaborators with a small sample 
of spiral (and the giant elliptical M87), but extends the idea to a 
much larger and diverse sample of shallow and steep cusp ellipti- 
cals and also allows for a direct parametric description of the 2D 
surface brightness profiles. 

We model the surface brightness (SB) profiles of 23 ellipti- 
cals in and around the Virgo cluster with a multi-component DW- 
function and summarize our observations as follows: 

1. Multi-component DW-function fits the SB profiles of el- 
lipticals with residuals consistently comparable to measurement er- 
rors over large dynamic ranges ~ 10^ in radii, and ~ 10** in SB, 
with a median sample rms of 0.032 mag arcsec^^. Nine galax- 
ies are well described with a 2-component (central and outer) DW 
model, while for fourteen galaxies a third component was required, 
and confirmed through an F-test. The third component acts as an in- 
termediate component between the central and outer components, 
except in NGC462 1 and NGC4434 where it appears to be embed- 
ded within the outer component. 

2. All steep-cusp and shallow-cusp galaxies reveal a central 
component that is in excess to an inward extrapolation of the outer 
component. Its shape parameter n is usually less than that of the 
outer component, which implies, as expected, that the central com- 
ponent is more concentrated than the outer component. Exceptions 
are NGC4459, NGC4387 and NGC4467 (sectionlOl. 

3. The central component of all massive shallow-cusp galax- 
ies have n < 1, while those of steep-cusp galaxies generally have 
n > 2, although there are some cases with 1 < n < 2 (section [83] l. 
This indicates that the central components of shallow-cusp galaxies 

a) are more concentrated than that of the steep-cusp galaxies 
even though the later are more denser in their central regions; and 

b) could signal the presence of disk-like systems; however, 
this must be verified spectroscopically. 

4. The central component of the shallow-cusp galaxies is far 
more luminous (and massive) and spatially more extended (large 
r_2 or rss) than that of the steep-cusp galaxies. Further, within 
a factor of two, all galaxies appear to host a similar fraction of 
total light in their central-l-intermediate components, with a weak 
indication that massive galaxies may be hosting a larger fraction. 
The last point is inconclusive due to the small number of massive 
ellipticals in our sample. 

5. In most of the shallow-cusp galaxies the outer compo- 
nent makes a comparable contribution to the density in the central 
regions with respect to that of the central component, while in the 
steep-cusp galaxies the central component is dominant. 

From the modelling point of view (section [83}, the shallow- 
cusp feature usually seen in massive galaxies is due to a combina- 
tion of a larger r_2, a low n < 1 implying a more concentrated 
component, and the non-negligible contribution to the density in 
the central regions from the n > 5 outer component. 

6. Galaxy formation models indicate that in massive ellipti- 
cals, the formation, evolution and subsequent coalesence of binary 
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SMBHs can remove nearly (2-3)Mbh of stellar material from the 
central regions, in a typical merger (refer to iii, section lll.2.3t . 
leading to mass deficits. Our observation that the central compo- 
nent in shallow-cusp galaxies are already massive, can be used to 
constrain galaxy formation models by accounting for the amount 
of mass that must have assembled in these galaxies prior to mass 
ejection by the binary SMBH's. 

It is possible that the current shapes of the central as well as 
the intermediate DW-components, within the central kiloparsec or 
so, have been influenced by the evolution of the binary SMBHs; 
and consequently some deficit could well have occured. However, 
the amount and sign of such deficits depend on estimating a break 
radius Rt, and assuming that the functional form of the density 
profile beyond Rt can be extrapolated all the way to the centre. 
This may not be meaningful and has led to large disagreements in 
the literature - both for a given galaxy as well as for averages over a 
sample. We note that core-scouring by SMBHs is unlikely to be the 
sole mechanism for producing some of the largest 'cores' (shallow- 
cusps), and other processes in galaxy formation and evolution are 
likely to have played their role as well in forming this feature. In 
such a case extrapolating the density profile inwards from Rt, to 
estimate the mass deficit, will yield misleading results. 

7. For 14 of the 23 galaxies we could describe the intrin- 
sic 3D luminosity density distribution fairly uniquely with a multi- 
component Einasto model (sectionfTOt. Since these galaxies span a 
wide range of luminosity — 24<AfvT< — 15 and come from both 
the steep- and shallow-cusp families, it is likely that other galaxies 
also have an intrinsic structure that can be modelled with a multi- 
component Einasto profile. 

8. Pure dark matter haloes are known to be well fit with the 
Einasto function. We have shown here that the same fitting function 
that describes the intrinsic density distribution of ACDM N-body 
haloes can also be used to model the intrinsic baryonic density of 
ellipticals, but for baryons a multi-component Einasto model is re- 
quired. There thus appears to be an universality in the functional 
form of the density profile of baryons (stars) and dark matter. 

9. The Einasto shape parameter n of dark matter haloes in N- 
body simulations is very similar to the n of the outer component of 
our massive ellipticals; both are 5 < n < 8. Further, PNe and strong 
lensing observations indicate that massive ellipticals are more dark 
matter dominated than less massive ellipticals. This indicates that 
our result - the outer component of the surface brightness profiles 
of massive galaxies has 5 < n < 8 - could imply i) a common 
feature of collisionless systems; and ii) that galaxies with such n 
for their outer (major) component are dark matter dominated. 

10. In section [TT] we have shown that the shapes of cen- 
tral and global profiles of shallow- and steep- cusp galaxies differ 
markedly and could be a result of differing formation pathways. 
In section [831 we have shown that the shallow-cusps can be typ- 
ically described as systems with - i) central DW-components of 
low n (< 1) and large scale radius r_2 or rsE, and ii) outer DW- 
component of large n (> 5) and central density comparable to that 
of the central DW-component. And in section [12] we have shown 
that large outer n systems are likely to be dark matter dominated. 

Hence, if dark matter has played a role in the formation and 
evolution of massive shallow-cusp galaxies, it will be instructive to 
explore its role, if any, in shaping the central regions and conse- 
quently in forming the shallow cusps. 

Finally we note that the galaxies modelled in this paper are all 
in and around the Virgo Cluster. Galaxy structure, formation and 



evolution is known to depend on the environment, and even though 
we suspect that our conclusion about the multi-component Einasto 
structure should be more widely applicable than Virgo, only more 
data can confirm that. 
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